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Introduction 

This  research  project  concerns  the  development  of  MRI  arterial  spin  tagging  to  non-invasively  measure 
breast  tissue  perfusion.  MR  dynamic  first-pass  contrast-enhanced  imaging  has  shown  that  malignant  and 
benign  breast  lesions  can  be  distinguished.  However,  it  may  have  limited  importance  in  clinical  breast 
diagnosis  due  to  significant  false-negatives  and  false-positives.  The  arterial  spin  tagging  technique  was 
developed  to  measure  tissue  perfusion  parameters  without  intravenous  contrast,  and  has  been  successfully 
demonstrated  in  brain  and  kidney.  The  specific  aims  are  to  (1)  refine  arterial  spin  tagging  pulse  sequences 
and  imaging  protocols,  (2)  develop  automated  data  analysis  software  for  measurement  of  breast  tissue 
parameters,  and  (3)  compare  the  technique  to  first-pass  contrast-enhanced  MRI  and  biopsy.  We  will  test 
the  hypothesis  that  arterial  spin  tagging  provides  accurate  and  precise  discrimination  between  normal 
tissue,  benign  and  malignant  lesions,  when  differences  in  perfusion  and  7j  exist.  Lesions  will  have  been 
previously  detected  by  clinically  accepted  diagnostic  imaging  procedures,  and  by  biopsy.  Statistical 
analysis  will  be  performed  to  access  the  correspondence  between  arterial  spin  tagging  and  biopsy,  and  to 
establish  the  relative  value  of  spin  tagging  compared  to  first-pass  contrast-enhanced  MRI.  We  hope  to 
establish  that,  relative  to  using  first-pass  contrast-enhanced  imaging,  false  positives  and  negatives  are 
reduced  using  arterial  spin  tagging  by  virtue  of  increased  image  signal-to-noise  ratio  (SNR),  higher  spatial 
resolution,  and  the  unique  ability  to  obtain  estimates  of  macromolecular  bound  fluid  fractions.  The  scope 
of  effort  on  the  project  is  mainly  limited  to  the  technical  aspects  of  development  of  a  new  methodology. 
However,  the  project  also  includes  a  rigorous  performance  comparison  with  the  current  gold-standard 
methodology.  The  technique  will  be  evaluated  in  sixty  patients,  with  roughly  equal  numbers  of  benign 
and  malignant  lesions. 
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Body 

Summary  of  Project  Accomplishments 

Report  for  2000 

David  Zhu  completed  his  dissertation  (attached  as  Appendix)  and  compiled  the  software  on  CD-ROM 
(enclosed),  a  poster  was  presented  at  the  BCRP  meeting  in  Atlanta  (attached  as  Appendix),  and  a 
publication  is  currently  being  written.  Comparison  of  the  arterial  spin  tagging  sequence  with  the  dynamic 
contrast  enhanced  first  pass  study,  intended  for  year  3,  has  not  yet  been  achieved.  Sixty  patients  should 
have  been  studied  at  this  point  in  the  project,  but  only  two  have  been  studied.  It  was  not  possible  to  find  a 
replacement  physician  for  the  original  co-investigator  on  the  project,  Rebecca  Zulim,  who  left  UCDMC 
November  1,  1998. 1  submitted  a  revised  project  schedule  and  budget  in  May  1999  based  upon  having  Dr. 
Philip  Schneider  as  a  collaborator,  and  these  were  approved,  but  they  were  not  implemented  due  to  the 
lack  of  patient  recruitment  into  the  study.  Although  Dr.  Schneider  agreed  to  serve  as  a  collaborator  for 
this  period,  patient  recruitment  was  ineffective.  Over  the  past  year,  exactly  $12,971  was  spent  from  grant 
funds,  for  support  of  Dr.  Zhu  during  the  last  three  months  (Oct  1999  -  Dec  1999)  of  his  employment  at 
UC  Davis,  and  3%  PI  salary  support.  No  salary  support  was  provided  to  the  co-I.  As  of  Sept  2000, 
$71,408  (27.2%  of  total  award)  remain  unexpended. 

From  the  1999  Annual  Report  . . . 

UCD  Medical  Center  Radiology  Department  upgraded  all  their  MRI  systems  from  Signa 
Genesis  5.x  to  Signa  LX  8.2.5,  resulting  in  several  months  during  which  our  pulse 
sequences  were  not  compatible  with  any  of  the  MRI  systems  to  which  we  had  access. 
Substantial  time  was  needed  to  convert  the  pulse  sequences  to  the  LX  8.2.5  platform 
from  the  Genesis  5.x  platform.  This  time  period  for  upgrading  the  UCD  Med  Center  MRI 
systems  was  not  planned  in  advance,  and  the  time  to  convert  our  research  software  was 
not  written  into  the  approved  Work  Statement.  The  graduate  student  on  this  project, 
David  Zhu,  and  I  succeeded  in  converting  the  arterial  spin  tagging  pulse  sequences, 
and  on  further  developing  new  arterial  spin  tagging  sequences  based  on  our  "odd- 
hybrid"  EPI  technique  [Buonocore  MH,  Zhu  DC.  High  spatial  resolution  EPI  using  an  odd 
number  of  interleaves.  Magnetic  Resonance  in  Medicine  41  (6):  1199-1205  (1999)]. 

Also,  advances  have  been  made  in  our  image  processing  software  for  the  project.  The 
major  deficiency  in  our  progress  is  with  the  third  technical  objective.  Comparison  of 
the  arterial  spin  tagging  sequence  with  the  dynamic  contrast  enhanced  first  pass 
study  has  not  yet  been  achieved.  Rebecca  Zulim,  the  co-I  on  the  project,  left  UCDMC 
November  1,  1998,  and  a  replacement  was  not  immediately  found.  The  budget  and 
schedule  was  officially  pushed  back  6  months,  into  the  3rd  year  which  had  a  light 
work  schedule  in  the  original  proposal.  The  project  now  continues  into  the  3rd  year, 
as  dictated  by  approved  revised  schedule  and  budget.  Dr.  PD  Schneider  of  the  UC  Davis 
Cancer  Center  is  the  co-investigator  on  the  project. 

From  the  1998  Annual  Report  ... 

There  are  two  substantial  contributions  of  the  work.  The  first  is  the  development  of 
robust  pulse  sequences  and  analysis  methods  to  derive  accurate  T1  estimates  with 
inherently  noisy  and  confounded  signal  data.  The  second  is  the  development  of  the 
software  program  (BrView) ,  to  allow  rapid  review,  quantitative  analysis,  and 
assessment  of  the  multitude  of  different  breast  images  and  timeseries  data  that  is 
obtained  in  each  patient  study. 
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Accomplishments  relating  to  Approved  Statement  of  Work 
Technical  objective  1 

Task  1;  Months  1-6:  Implementation  and  testing  of  magnetization  transfer  pulses  for  both 
arterial  tagging  and  first-pass  contrast  enhanced  sequences. 

Report  for  2000 

Magnetization  transfer  (MT)  was  fully  tested  with  the  EPI-based  spin  tagging  sequence. 

From  1999  Annual  Report  ... 

Magnetization  transfer  (MT)  has  been  implemented  with  the  EPI-based  spin  tagging 
sequence  (see  Figure  1) .  The  MT  radio frequency  (RF)  pulse  selectively  excites  the 
protons  of  the  macromolecules  and  its  hydration  layer  by  a  900  flip  angle.  Two 
gradient  spoiling  pulses  (in  the  x-  and  y-  directions)  are  applied  immediately  after 
the  RF  pulse  to  spoil  the  transverse  magnetization.  The  usual  RF  excitation  pulse  for 
data  acquisition  is  applied  immediately  after  these  spoiling  pulses.  Due  to  hardware 
and  or  compiler  problems  that  are  currently  being  investigated  by  GE,  this  pulse 
sequence  has  been  implemented  only  on  pulse  sequence  simulation  software  (EPIC) ,  not 
yet  on  the  MR I  system.  The  hardware  problems  are  somehow  preventing  this  sequence 
from  running  successfully  on  the  MRI  system.  Testing  of  MT  using  the  fast  SPGR-based 
sequence  revealed  that  it  would  not  be  effective.  Because  the  magnetization  transfer 
RF  pulse  is  long  duration  (16  ms) ,  the  "small  TR"  (optr)  period  of  the  fast  SPGR 
sequence  was  effectively  doubled  in  length,  thus  sacrificing  short  scan  time  and 
distorting  the  signal  upon  which  the  perfusion  measurement  is  based. 

From  1998  Annual  Report ... 

Magnetization  transfer  pulses  have  not  been  implemented  yet.  Mainly,  this  delay  was 
due  to  the  fact  that  the  MRI  system  at  UC  Davis  Medical  Center  is  scheduled  to  be 
upgraded  to  the  new  Horizon  LX  system  soon,  and  we  opted  to  defer  this  pulse  sequence 
development  until  we  were  using  this  new  LX  software  platform.  There  have  been 
multiple  delays,  but  this  upgrade  is  scheduled  to  occur  by  January  1999.  The  existing 
pulse  sequence  will  be  converted  from  the  5.4  OS  platform  to  the  LX2  platform  (a 
substantial  change) ,  and  the  magnetization  transfer  pulses  will  be  implemented. 

Because  of  concerns  regarding  the  ultimate  utility  of  the  magnetization  transfer 
technique,  we  have  decided  that  implementation  and  testing  of  arterial  spin- tagging 
sequence  based  on  echo  planar  imaging  (EPI)  data  acquisition  is  a  higher  priority. 

EPI  based  spin-tagging  technique  can  theoretically  improve  the  accuracy  of  the  2\ 
estimation,  which  we  have  studied  and  worked  with  extensively  (see  below) .  The  major 
problem  with  magnetization  transfer  is  that  it  represents  a  perturbation  on  an 
already  small  signal.  Therefore,  whether  we  will  observe  a  magnetization  transfer 
effect  is  questionable.  Nevertheless,  it  will  be  implemented  as  part  of  the  project. 

Task  2:  Months  1-9:  Implementation  and  testing  of  interleaved  high-resolution  imaging 
technique,  for  both  arterial  tagging  and  first-pass  contrast  enhanced  sequences. 

Report  for  2000 

Conversion  of  the  Fast  SPGR  Based  Arterial  Spin  Tagging  Sequence 
The  authors  completed  the  conversion  and  testing  of  the  fast  SPGR-based  arterial  spin  tagging  pulse 
sequence,  from  the  Genesis  5.4  platform  to  the  LX  8.2.5  platform. 

Development  of  an  EPI-Based  Arterial  Spin  Tagging  Sequence 


6 


In  October  1999,  the  authors  formulated  an  alternative  approach  to  eliminate  the  requirement  of  polarity 
alternation,  which  was  successfully  implemented.  The  following  steps  are  performed  at  each  interleaf:  (1) 
non-selective  RF  inversion  pulse,  (2)  20  EPI  acquisitions,  each  of  which  is  preceded  by  a  10°  RF  pulse 
(Figure  2).  These  20  acquisitions  are  played  out  with  a  repetition  time  ( optr  in  Figure  2)  of  100  ms.  Each 
acquisition  defines  one  point  on  a  T\  recovery  curve.  The  total  time  for  each  interleaf  (TR  in  Figure  5.10) 
is  2  seconds.  Data  for  the  other  interleaves  is  collected  in  the  same  fashion  in  subsequent  TR  intervals  (see 
Dissertation).  The  completed  set  of  images,  showing  T1  recovery  of  the  magnetization,  is  acquired  in  two 
seconds  times  the  number  of  interleaves. 

From  1999  Annual  Report . . . 

Summary 

Since  the  completion  of  the  upgrade  of  the  UC  Davis  Imaging  Center  MR I  system,  and 
since  it's  designation  as  the  Research  MR  I  system  for  UC  Davis,  we  have  thoroughly 
studied  the  EPIC  LX  8.2.5  programming  language  for  pulse  sequence  development,  and 
the  operation  of  this  new  MRI  system.  The  transition  to  the  new  research  MRI  required 
a  complete  rewrite  of  all  pulse  sequence  software  that  we  had  developed  for  this 
project.  In  addition,  we  revised  the  pulse  sequences  to  utilize  the  high  performance 
gradients  (40  mT/m  peak,  150  mT/m/ms  rise)  of  the  system.  We  believe  these  revised 
sequences  will  provide  better  arterial  spin  tagged  data  for  perfusion  measurement. 

All  of  the  other  sequences  in  our  breast  imaging  protocol  were  set  up  for  the  new  LX 
8.2.5  platform  as  well. 

Conversion  of  the  Fast  SPGR  Based  Arterial  Spin  Tagging  Sequence 
The  authors  converted  the  fast  SPGR-based  arterial  spin  tagging  pulse  sequence  from 
the  Genesis  5.4  platform  to  the  LX  8.2.5  platform.  The  converted  sequence  was  also 
improved  to  achieve  higher  resolution  with  a  shorter  overall  scan  time.  The  new  MR 
system  has  higher  peak  gradient  amplitude  and  slew  rate.  Fast  data  acquisition,  which 
requires  the  higher  read-out  gradient  amplitude  (4.0  g/cm)  could  be  specified,  and 
consequently  the  data  acquisition  period  (optr)  was  shorten  from  16.62  ms  to  12.5  ms. 
Since  optr  was  shortened,  a  larger  matrix  in  the  phase  encode  direction  could  be 
specified.  The  resolution  is  now  256  x  256  matrix  using  a  one-half  phase  field  of 
view  acquisition.  Previously  we  used  a  256  x  240  matrix  with  one-half  phase  field  of 
view  acquisition.  The  total  time  for  the  image  acquisition  (the  TR  period) ,  including 
the  global  inversion  RF  pulse,  has  been  decreased  from  2.7  seconds  to  2.22  seconds. 
The  performance  of  this  new  sequence  has  been  confirmed  using  phantoms  and  one  test 
subject . 


Development  of  an  EPI -Based  Arterial  Spin  Tagging  Sequence 
The  EPI-based  arterial  spin  tagging  sequence  was  rewritten  and  verified  on  the 
research  MRI  system.  The  data  acquisition  scheme  is  based  on  our  recently  published 
high-resolution  odd-number  interleaf  EPI  sequence  [Buonocore  MH,  Zhu  DC.  High  spatial 
resolution  EPI  using  an  odd  number  of  interleaves.  Magnetic  Resonance  in  Medicine  41 
(6):  1199-1205  (1999)].  This  sequence  was  first  developed  and  verified  on  the  Signa 
Advantage  1 . 5  T  GE  MR  system.  Unfortunately,  under  the  new  LX  system,  the  polarity 
alternation  scheme  for  each  successive  interleaf  we  had  implemented  caused  the  system 
to  hang,  and  substantial  time  was  expended  understanding  the  problem. 

From  1998  Annual  Report  ... 

An  interleaved  arterial  spin-tagging  technique,  that  allows  high  resolution  (e.g.  256 
x  256  or  more)  has  been  implemented  and  tested.  We  also  implemented  a  rectangular 
data  acquisition  technique,  which  acquires  data  sets  that  have  different  numbers  of 
points  along  the  frequency  and  phase  encode  direction.  The  256  x  240  interleaved 
acquisition  provides  the  best  spatial  versus  temporal  resolution  tradeoff,  and  has 
been  used  for  all  of  our  recent  studies. 


7 


Task  3;  Months  1-9:  Analysis  of  spin  tagging  sequence  to  understand  causes  of  existing  baseline 
offsets,  effects  of  inversion  slice  transition  profiles,  and  effects  of  RF  flip  angle  profile  on  the 
measurements,  with  implementation  and  testing  pulse  sequence  modifications  to  minimize  these 
imperfections. 

Report  for  2000 

The  EPI  arterial  spin  tagging  sequence  was  not  formally  evaluated. 

From  1999  Annual  Report  . . . 

Regarding  the  spin  tagging  sequence  utilizing  fast  SPGR,  these  issues  were  treated  in 
the  1998  report.  The  EPI  arterial  spin  tagging  sequence  has  not  been  evaluated,  but 
the  same  solutions  for  the  inversion  slice  transition  profiles,  and  effects  of  RF 
flip  angle  profile,  will  apply. 

From  1998  Annual  Report  . . . 

Matching  of  inversion  and  excitation  slice  profiles  was  greatly  improved  by  using 
customized  RF  pulses  designed  with  the  Shinnar-Le  Roux  (SLR)  algorithm.  The  new  data 
processing  technique  (based  on  a  semi-log  linear  regression  of  inversion  time  (TI) 
dependent  signal)  significantly  reduced  the  previously  reported  problem  in  measuring 
the  longitudinal  relaxation  time  (Tl)  caused  by  baseline  offsets,  and  by  the 
uncertainty  in  the  effective  inversion  time.  The  optimal  RF  flip  angle  for  spin 
excitation  during  spatial  encoding  was  found  to  be  10  degrees,  based  on  experiments 
over  a  range  of  flip  angles. 

Technical  objective  2 

Task  4:  Months  3-15:  Software  for  automatic  estimation  and  error  analysis  of  perfusion,  tissue 
water  longitudinal  relaxation  time,  and  extracellular  fluid  volume  fraction  from  mathematical 
models  and  user-defined  ROIsfrom  spin  tagging  timeseries.  Implementation  of  pharmacokinetic 
model  calculations  and  error  analysis  for  first-pass  contrast  enhanced  imaging. 

Report  for  2000 

No  new  software  was  developed. 

From  1999  Annual  Report . . . 

From  October  1998  through  April  1999  (during  which  time  there  was  no  MRI  system  for 
us  to  develop  on,  or  do  patient  studies) ,  image  processing  software  development 
continued.  Improvements  in  the  BreastView  program,  development  of  "dispAlls"  program 
for  analysis  of  general  MRI  images,  and  implementation  of  statistical  analysis  of 
breast  lesions  based  on  Bayesian  approach,  was  completed.  Statistical  analysis  of 
available  breast  lesion  data  continued  (See  below) .  The  most  important  new  function 
of  BreastView  is  that  for  analysis  of  the  first-pass  contrast  enhanced  images.  The 
program  creates  an  image  based  upon  the  time  derivative  of  the  rate  of  rise  of  the 
signal,  so  that  the  contrast  enhancement  effect  can  be  better  visualized  (See  Figures 
4  and  5) . 


From  1998  Annual  Report  ... 

A  software  program  (BrView)  has  been  written  in  C,  X  Window  System,  and  Motif,  and 
implemented  on  an  SGI  02  computer  (paid  for  by  grant  funds)  for  analyzing  and 
visualizing  the  MR  images  and  timeseries.  Much  effort  has  been  expended  on  developing 
the  capability  to  easily  review  all  of  the  images  taken  in  each  study,  and  cross- 
reference  pixel  locations.  Much  more  effort  than  anticipated  was  required  for 
development  of  a  robust  technique  for  identifying  so-called  "suspicious  regions" 
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based  on  the  T1  and  perfusion  measurements.  The  measurement  of  Tl,  upon  which  the 
perfusion  measurement  entirely  depends,  is  based  on  a  semi-log  linear  regression 
technique  developed  by  the  investigators.  The  Tl,  perfusion  (f/A.)  and  standard  errors 
of  these  quantities  are  estimated  automatically  using  this  robust  technique.  The  so- 
called  feature  images  display  these  quantities  in  color  and  grayscale.  These 
quantities  are  used  to  calculate  a  "suspicion  index"  for  each  pixel,  and  thereby 
identify  regions  of  breast  tissue  suspicious  for  malignancy. 

The  implementation  of  the  pharmacokinetic  model  calculations  and  error  analysis  for 
first-pass  contrast  enhanced  imaging  has  not  been  completed.  Thus  far,  the  software 
allows  the  user  to  click  on  the  reference  image  of  the  dynamic  study,  to  display  the 
time  profile  of  the  signal  at  that  pixel.  The  implementation  of  the  pharmacokinetic 
model  is  currently  a  major  focus  of  the  current  effort. 

Task  5:  Months  3-15:  Software  for  registration  (including  implementation  and  testing  of  motion 
correction  and  physiological  noise  reduction  algorithms)  and  overlay  of  images  from  high 
resolution  T /  weighted,  spin  tagging,  and  contrast  enhanced  studies. 

Report  for  2000 

BreastView 

BreastView  now  performs  a  simple  analysis  of  first -pass  contrast  enhanced  images.  The  program  creates  a 
contrast  enhancement-weighted  image  based  upon  the  maximum  time  derivative  of  the  rate  of  rise  of  the 
signal,  so  that  the  contrast  effect  can  be  better  visualized. 

The  dispAlls  program 

There  were  no  further  enhancements  of  this  program. 

From  the  1999  Annual  Report . . . 

First  described  in  the  1998  Annual  Report,  the  BreastView  program  has  been 
continually  upgraded  and  enhanced  with  the  latest  algorithms  for  processing  and 
analysis  of  breast  imaging  studies,  including  those  for  spin-tagged  and  dynamic 
contrast  enhanced  images.  The  specific  developments  this  year  of  BreastView,  and  of 
another  program,  dispAlls,  are  described. 

BreastView 

Between  October  1998  and  April  1999,  the  BreastView  program  was  made  "clinician- 
friendly"  .  The  graphical  user  interface  now  guides  the  user  through  the  entire 
visualization  and  analysis  process  for  a  breast  imaging  study.  A  "Help"  button, 
linked  to  extensive  set  of  manual  pages,  now  appears  on  the  first  page  of  the 
interface.  Clicking  the  Help  button  brings  up  an  html  file  with  necessary  operating 
instructions.  The  help  files  were  written  for  clinicians.  The  clinician  can  now  begin 
using  the  program  without  first  studying  the  written  instructions.  If  an 
inappropriate  option  is  selected,  the  program  now  provides  help  to  the  user  to  select 
a  more  appropriate  option  or  operation.  Several  new  features  and  functions  have  been 
implemented.  The  most  important  new  function  of  BreastView  is  the  analysis  of  the 
first-pass  contrast  enhanced  images.  The  program  creates  a  contrast  enhancement- 
weighted  image  based  upon  the  time  derivative  of  the  rate  of  rise  of  the  signal,  so 
that  the  contrast  effect  can  be  better  visualized. 

The  dispAlls  program 

This  program  has  been  written  for  quickly  viewing  and  performing  ROI  analysis  on  one 
or  a  series  of  images.  The  program  is  particularly  useful  for  quickly  reviewing  the 
timeseries  of  signal  changes  in  an  ROI  from  the  dynamic  contrast  enhanced  studies. 

The  program  displays  one  image,  or  a  sequence  of  images,  with  or  without  an  image 
header.  It  performs  ROI  analysis  and  plots  the  calculated  values. 
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From  1998  Annual  Report  . . . 

Using  the  BrView  software  program,  image  pixels  containing  high  values  of  perfusion 
( f/X ),  and  moderate  to  high  values  of  Tlt  and  low  standard  errors,  are  identified  as 
suspicious,  and  assigned  a  "suspicion  index"  based  on  finding  similar  abnormal  values 
at  spatially  adjacent  pixels.  Suspicious  pixels  can  be  overlaid  on  the  high- 
resolution  Ti  and  T2  clinical  images  by  color  mapping,  and  simultaneously  presented 
with  the  first-pass  contrast  enhanced  timeseries  at  these  pixels.  This  provides  the 
user  with  a  comprehensive  anatomical  and  functional  view  of  the  suspicious  regions, 
and  facilitates  making  a  decision  regarding  the  malignant  nature  of  the  lesion. 
Finally,  a  motion  artifact  estimation  and  correction  algorithm  has  been  developed  and 
implemented  in  a  separate  software  program. 

Technical  objective  3 

Task  6:  Months  9-24:  60  patients  with  malignant  and  benign  breast  lesions  will  be  imaged  using 
T/  weighted,  first-pass  contrast  enhanced,  and  arterial  spin  tagging  MRI  pulse  sequences. 

Report  for  2000 

No  additional  subjects  have  been  done. 

From  the  1999  Annual  Report . . . 

Unfortunately,  no  additional  subjects  have  been  done.  This  is  due  to  the  loss  of 
Rebecca  Zulim,  MD,  as  the  co-investigator,  and  the  transition  to  a  dedicated  research 
MRI  system  that  operated  at  LX  8.2.5,  as  described  above.  Through  Sept  1998,  eighteen 
subjects,  including  11  patients  referred  by  their  physicians  and  seven  self-referred 
volunteers,  had  participated  in  the  breast  studies.  Three  patients  and  two  volunteers 
were  eliminated  from  the  analysis  because  the  studies  were  not  done  according  to 
protocol,  due  to  excess  motion,  or  the  abnormal  mass  being  too  close  to  the  chest 
wall  to  be  measured  reliably,  or  no  patient  biopsy  to  confirm  the  MR  study  result. 
There  are  total  of  13  useful  cases  left  for  analysis.  For  all  these  cases,  the 
suspicious  pixels  were  identified  based  on  the  following  criteria:  Tin:  >  0.5  sec 
f/DD.>  0.1  sec-1,  STD  of  f/D:  <  0.1  sec-1,  and  the  suspicion  level  threshold  was  set 
at  20.2%.  Analysis  has  been  done  using  the  program  BreastView.  Table  8.5  shows  the 
result  of  these  13  cases. 

From  1998  Annual  Report  ... 

We  have  thus  far  recruited  two  subject  studies  for  both  first-pass  contrast  enhanced 
and  arterial  spin  tagging  MRI  pulse  sequences.  Now  that  the  software  development  is 
substantially  completed,  we  are  trying  to  improve  the  rate  at  which  patients  are 
recruited  into  the  study.  It  is  likely  that  the  patient  studies  will  extend  into  the 
third  year  of  the  proposal.  We  are  currently  discussing  with  physicians  at  a  local 
Kaiser  Hospital  for  recruitment  of  additional  patients.  Appropriate  documentation 
will  be  submitted  for  approvals  upon  reaching  a  firm  agreement. 


Task  7  (listed  as  Task  6  in  the  original  grant  application):  Months  12-36:  Automated 
pharmacokinetic  analysis,  blinded  image  reading,  and  statistical  comparison  of  arterial  spin 
tagging,  contrast  enhanced  MRI,  and  biopsy  results. 

Report  for  2000 

No  new  software  was  developed. 
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From  the  1999  Annual  Report  . . . 

We  developed  an  analysis  program  for  the  detection  of  the  features  of  the  dynamic 
contrast  enhanced  study,  and  have  developed  a  statistical  approach  to  the 
determination  of  the  malignant  status  of  the  lesion. 

Bayesian  Statistical  Analysis 

The  identification  of  a  suspicious  pixel  based  upon  the  index  threshold,  that  was 
developed  last  year,  is  a  simple  approach  to  cancer  detection.  When  a  large  number  of 
studies  become  available,  a  Bayesian  decision  technique  [53]  will  be  applied.  This 
technique  will  be  used  to  provide  ROC  (receiver  operating  characteristic)  curves  that 
are  commonly  used  for  evaluation  of  diagnostic  algorithms.  To  summarize  the  Bayesian 
technique,  let  vector  space  x  =  (Tl,  f/X,  STD  of  f/X,  etc.)  for  a  ROI  (region  of 
interest)  in  a  breast. 

From  the  1998  Annual  Report  . . . 

Automated  pharmacokinetic  analysis,  blinded  image  reading,  and  statistical 
comparisons,  have  not  been  started. 

Key  Research  Accomplishments 

Report  for  2000 

A.  EPI-based  arterial  spin  tagging  sequence,  with  MT  pulse  option,  for  Signa  Horizon  LX  CV/i  system. 
From  the  1999  Annual  Report  . . . 

A.  Pulse  sequences  that  are  compatible  with  a  Signa  Horizon  LX  EchoSpeed  CV/i  system, 
running  LX  8.2.5  system  software: 

1.  Fast  SPGR-based  arterial  spin  tagging  sequence 

2 .  Odd-number  interleaved  EPI  sequence 

3.  Interleaved  EPI-based  arterial  spin  tagging  sequence 

B.  Imaging  processing  software 

1 .  BreastView 

2.  DispAlls 

3 .  Bayesian  analysis  of  arterial  spin  tagged  image  data 

From  the  1998  Annual  Report  ... 

A.  Fast  SPGR-based  arterial  spin  tagging  sequence  for  Signa  Advantage  MRI  system, 
running  5.7  system  software: 

B.  Imaging  processing  software  (BrView) 

Reportable  Outcomes 

Report  for  2000 

A.  Buonocore  MH,  Zhu  DC.  Magnetic  resonance  arterial  spin  tagging  for  non-invasive 

pharmacokinetic  analysis  of  breast  cancer.  Poster  Presentation  and  Abstract.  Proceedings  of  the 
Era  of  Hope,  Dept,  of  Defense  Breast  Cancer  Research  Program  Meeting,  June  8-11,  2000,  Hilton 
Atlanta  and  Towers,  Atlanta,  GA,  Vol.  1,  page  177. 


11 


B.  Buonocore  MH,  Zhu  DC,  Zulim  RA.  Magnetic  resonance  arterial  spin  tagging  for  non-invasive 
pharmacokinetic  analysis  of  breast  cancer.  Poster  Presentation  and  Abstract.  Proceedings  of  the 
UC  Davis  Cancer  Center  Symposium,  Oct  6-7,  2000.  UC  Davis  Cancer  Center,  Sacramento,  CA. 

C.  David  Zhu.  PhD  Dissertation.  Magnetic  resonance  pulse  sequences  and  analytical  techniques  for 
breast  cancer  detection  and  cardiovascular  flow.  UC  Davis,  Biomedical  Engineering.  Dec  1999. 

D.  David  Zhu.  PhD  in  Biomedical  Engineering,  Dec  1999. 

E.  David  Zhu,  PhD.  Employment  as  GE  Applications  Engineer,  beginning  Jan  2000. 

From  the  1999  Annual  Report  . . . 

A.  Buonocore  MH,  Zhu  DC.  High  spatial  resolution  EPI  using  an  odd  number  of 
interleaves.  Magnetic  Resonance  in  Medicine  41  (6):  1199-1205  (1999). 

B.  Buonocore  MH,  Zhu  DC,  Zulim  RA.  Analysis  software  for  breast  imaging  studies. 
Proceedings  of  the  International  Society  for  Magnetic  Resonance  in  Medicine, 
7th  Annual  Meeting  and  Exhibition,  3:  2172  (1999). 

C.  Pulse  sequences  for  Version  8.2.5  (MH  Buonocore  and  D  Zhu) : 

Arterial  spin  tagging  sequence  with  Fast  SPGR  Acquisition 
Odd-number  interleaf  EPI 

Arterial  spin  tagging  sequence  with  Odd-number  interleaf  EPI  acquisition 
Arterial  Spin  Tagging  sequence  with  Even-number  hybrid  EPI  acquisition 

D.  Breast  Imaging  Protocols  (MH  Buonocore  and  D  Zhu) 

Anatomical,  spin  tagged,  and  dynamic  c-e  breast  imaging  protocol 

E.  SGI  Graphics  Software  for  Breast  image  analysis  (MH  Buonocore  and  D  Zhu) : 
BreastView  (display  program  for  breast  imaging  protocol) 

DispAll  (display  program  for  general  MR I  images) 


From  the  1998  Annual  Report  . . . 

A.  Buonocore  MH,  Zhu  DC.  Odd-number  hybrid  EPI.  Proceeding  of  the  International 
Society  for  Magnetic  Resonance  in  Medicine,  6th  Annual  Meeting  and  Exhibition, 
on  CD-ROM,  p.  1967  (1998). 

B.  Buonocore  MH,  Zhu  D,  Pellot-Barakat  C.  Measurement  of  breast  tissue  perfusion 
using  arterial  spin  tagging,  Proceedings  of  the  International  Society  for 
Magnetic  Resonance  in  Medicine  5th  Annual  Meeting  and  Exhibition,  1:  311 
(1997)  . 

C.  Buonocore  MH,  Zhu  DC,  Pellot-Barakat  C,  Zulim  RA.  Noninvasive  measurement  of 
blood  flow  through  breast  tumors.  Book  of  Abstracts,  California  Breast  Cancer 
Research  Symposium.  Sacramento  Convention  Center,  Sacramento  CA,  Sept.  16, 
1997. 

D.  Buonocore  MH,  Zhu  DL,  Pellot-Barakat  C.,  Zulim  RA.  Non-invasive  measurement  of 
breast  tissue  perfusion  using  arterial  spin  tagging.  Radiology,  November  1997, 
205  (P) :  162. 

E.  Sequence  and  Protocols  for  arterial  spin  tagging  for  Version  5.7  Signa 
Advantage  MR I 

F.  Software  for  display  and  analysis  of  breast  images,  written  in  C,  X  and  Motif. 
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Conclusions 


Report  for  2000 

Technical  accomplishments  over  the  past  three  years  have  been  documented  in  a  PhD  Dissertation,  a  CD 
ROM  of  research  software,  and  poster  presentations.  These  will  be  used  to  accomplish  the  third  Specific 
Aim  of  the  project,  to  non-invasively  measure  perfusion  parameters  in  breast  masses  prior  to  biopsy. 

These  measurements  may  distinguish  malignant  from  benign  breast  lesions. 

From  the  1999  Annual  Report . . . 

We  have  successfully  converted  and  enhanced  the  arterial  spin  tagging  pulse  sequences 
to  work  on  a  new  research  MRI  system,  and  have  made  several  enhancements. 
Understanding  of  the  8.2.5  pulse  sequence  language  on  the  new  system  represents  a 
major  advance  of  the  project.  These  pulse  sequences  will  be  used  in  all  future  work, 
and  hopefully  there  will  not  be  a  substantial  MRI  system  upgrade  such  as  we 
experienced  for  a  few  years  at  least.  Image  display  and  statistical  analysis  software 
has  been  written  for  the  spin  tagging  and  contrast  enhanced  studies,  and  to  allow 
analysis  to  be  done  by  a  busy  clinician  without  computer  expertise. 

From  the  1998  Annual  Report  . . . 

In  this  year  of  the  project,  we  have  developed  the  necessary  pulse  sequences  for 
anatomical  and  functional  imaging  of  breast  tissue,  and  the  necessary  software  for 
the  analysis,  display  and  interpretation  of  this  data.  These  developments  will  allow 
us  to  evaluate  complete  breast  imaging  studies  that  are  composed  of  anatomical  scans, 
dynamic  first  pass  contrast  enhanced  scans,  and  arterial  spin  tagging  scans.  The 
integrated  nature  of  the  software  allows  easy  and  consistent  interpretation  of  both 
spatial  dependencies  as  well  as  temporal  dependencies  of  signal  changes  indicative  of 
disease.  The  in-vitro  and  in-vivo  studies  done  thus  far  indicate  that  first  and 
foremost,  arterial  spin  tagging  is  a  sensitive  and  reliable  method  for  measuring  Tl 
and  parameters  related  to  tissue  perfusion.  Second,  they  indicate  that  arterial  spin 
tagging  is  a  viable  alternative  to  dynamic  first-pass  contrast-enhanced  imaging.  The 
fact  that  arterial  spin  tagging  can  be  easily  added  to  any  standard  breast  imaging 
protocol,  without  requiring  special  nursing  or  MR  technologist  expertise,  means  that 
it  is  especially  attractive  for  routine  clinical  use.  It  could  be  used  to  supplement 
contrast-enhanced  studies  in  the  initial  study  of  the  patient,  and  used  exclusively 
in  follow-up  studies.  Since  10/1/98,  six  subjects  have  been  studied  without  the 
dynamic  first-pass  contrast-enhanced  scan,  and  two  have  been  studied  with  this  scan. 
Comparison  of  contrast-enhanced  images  versus  arterial  spin-tagged  images  has  not  yet 
been  made.  In  particular,  no  statistical  comparison  of  the  dynamic  study  results 
(considered  to  be  the  gold  standard  for  non-invasive  assessment  of  malignant  disease) 
with  the  arterial  spin  tagging  results  has  been  performed.  However,  there  is 
considerable  spatial  overlap  of  suspicious  regions  obtained  by  these  methods  noted  on 
visual  inspection.  Collecting  data  on  human  subjects,  and  statistical  confirmation  of 
this  agreement,  will  be  major  focuses  of  next  year's  work. 
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Abstract 


This  dissertation  involves  the  development  of  magnetic  resonance  pulse 
sequences  and  analytical  techniques  for  two  projects:  breast  cancer  detection  and 
cardiovascular  flow  study. 

In  the  first  project,  arterial  spin  tagging  techniques  and  analytical  methods  are 
developed  and  implemented  to  measure  spin-lattice  relaxation  time  7j,  relative  blood 
perfusion,  and  magnetization  transfer  in  breast  tissue.  Studies  have  shown  that  water 
content  in  malignant  tumors  is  higher  than  in  normal  tissues.  Studies  have  also  shown 
that  perfusion  in  malignant  tumors  is  higher  than  in  both  normal  tissues  and  benign 
lesions.  Compared  to  the  current  widely  used  dynamic  contrast  enhancement 
techniques  for  breast  cancer  detection,  the  arterial  spin  tagging  techniques  may 
provide  more  sensitive  and  specific  assessment  of  tissue  parameters  at  higher  spatial 
resolutions  by  virtue  of  their  ability  to  improve  signal  to  noise  ratio  (SNR)  over 
reasonable  scan  times. 

>  >w 

For  the  second  project,  a  three-dimensional  (3D)  cine  phase  contrast  pulse 
sequence  (also  called  a  4D  sequence  in  this  dissertation)  is  developed  and 
implemented.  The  sequence  development  is  based  on  a  two-dimensional  (2D)  cine 
phase  contrast  technique  and  a  3D  non-cine  phase  contrast  technique  that  have  been 
developed  and  used  in  clinical  studies.  In  the  2D  technique,  a  sequence  of  images  is 
acquired  to  capture  the  motion  of  blood  flow  throughout  a  cardiac  cycle,  but  they  are 
limited  to  one  thin  slice  of  the  body.  Images  must  be  acquired  at  multiple  slice 
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locations  in  order  to  show  the  volume  change  of  flow  patterns.  The  3D  technique 
allows  the  acquisition  of  a  slab  of  images,  but  does  not  provide  the  dynamic 
information  of  cardiac  flow.  The  4D  technique  offers  the  convenience  of  volume 
imaging  of  the  3D  technique,  and  the  dynamic  cardiac  flow  information  of  the  2D 
technique.  The  flow  patterns  based  on  the  processed  images  can  be  used  to  assess  the 
cardiac  condition  of  the  subject. 
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Parti 


Arterial  Spin  Tagging  Techniques 

for  Breast  Cancer  Detection 

Chapter  3 
Introduction 

Breast  cancer  is  the  leading  cause  of  death  from  cancer  among  women  in  the 
western  world.  In  the  U.S.,  it  accounts  for  approximately  26%  of  all  cancer  cases  and 
1 8%  of  all  cancer  deaths  in  women.  Early  detection  and  treatment  have  played 
important  roles  in  improving  the  chance  of  survival  in  women  with  breast  cancer 

[18,19,20],  However,  early  detection  of  cancer  has  not  been  easy,  which  is  often 

1 

misled  by  far  more  frequent  incidences  of  benign  disorders.  The  ratio  of  benign  to 
malignant  lesions  is  greater  than  25:1  [18].  Few  benign  disorders  are  life-threatening. 
X-ray  mammography  is  currently  Used  as  the  primary  screening  tool  for  breast  cancer. 
If  suspicion  regions  were  found  with  X-ray  mammography,  a  patient  would  be  sent  to 
the  pathology  department  for  a  biopsy,  which  serves  as  the  final  measurement  of  the 
malignancy  of  the  suspicious  regions.  Although  it  is  crucial  to  detect  all  the  malignant 
cases,  biopsies  are  expensive  and  invasive.  It  is  economical  and  desirable  to  reduce  the 
number  of  benign  cases  sent  for  biopsy.  The  sensitivity  of  X-ray  mammography  is 
greater  than  90%.  But  the  percentage  of  suspicious  X-ray  findings  proven  to  be 
cancerous  (the  predictive  value  of  X-ray)  is  only  10-54%  [18].  Despite  its  low 
predictive  value  as  well  as  the  discomfort  of  breast  compression  and  potential 
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radiation  damage,  X-ray  mammography  is  convenient  and  economical.  The  cost  of  X- 
ray  mammography  is  about  $55  as  compared  to  about  $150  for  breast  ultrasound,  $150 
for  fine-needle  aspiration  of  palpable  masses,  and  $850  for  core  biopsy  of  occult 
lesions  [21].  Thus,  X-ray  mammography  will  continue  to  serve  as  an  important  role  in 
breast  screening. 

Due  to  the  drawbacks  of  X-ray  mammography,  researchers  have  studied  other 
imaging  modalities,  such  as  ultrasonography,  thermography,  computed  tomography 
(CT)  and  magnetic  resonance  imaging  (MRI)  [18, 22, 23].  Ultrasonography  is  an 
important  complementary  breast  imaging  tool,  which  can  make  the  possible 
differentiation  between  cystic  and  solid  masses  [20].  But  its  ability  in  detecting  tumors 
less  than  1  cm  and  non-palpable  malignancies  is  weak,  with  a  detection  accuracy  of 
only  about  30%  and  8%  respectively.  In  addition,  the  sensitivity  is  only  about  64% 
[18].  Theromography  is  a  modality  based  on  the  skin  thermal  image  after  the 

application  of  infra-red  radiation.  Although  this  method  has  been  postulated  as  a  risk 

* 

indicator  for  subsequent  development  of  breast  cancer,  it  has  not  been  proven  to  be 
useful  for  breast  cancer  detection  [20].  CT  has  appeared  to  be  more  sensitive  in 
detecting  breast  malignancies  than  X-ray  mammography  [18].  But  CT  of  the  breast 
requires  intravenous  injection  of  contrast  media  to  achieve  a  good  result.  Besides  this, 
the  cost  of  contrast  media,  in  addition  to  the  small  number  of  patients  examined  per 
hour  and  the  potential  radiation  damage,  have  limited  its  application  as  a  screening 
tool.  It  is  used  only  occasionally  as  a  diagnostic  test  [20],  MRI  is  new  in  breast 
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imaging.  But  its  great  sensitivity  to  soft  tissue  imaging  and  its  non-invasiveness  have 
provided  a  great  potential  for  breast  cancer  detection. 

In  MRI,  the  most  popular  technique  for  breast  cancer  detection  is  MR  dynamic 
first-pass  contrast  enhancement  (commonly  named  Magnetic  Resonance 
Mammography).  MR  dynamic  first-pass  contrast  enhancement  has  shown  its 
capability  in  distinguishing  between  malignant  and  benign  breast  lesions  [18]. 
However,  it  may  have  limited  importance  in  clinical  breast  diagnosis  due  to  significant 
false-negatives  and  false-positives  resulting  from  infectious  lesions,  proliferative 
dysplasia,  myxoid  fibroadenoma,  and  other  benign  breast  lesions.  Additionally, 
dynamic  contrast  enhancement  cannot  be  used  in  screening,  even  within  defined  high- 
risk  sub-populations  (e.g.  patients  with  dense  breasts),  due  to  the  high  cost  of  contrast 

I 

material  and  its  administration. 

An  arterial  spin  tagging  technique  J24]  has  been  developed  for  measuring 
tissue  parameters  without  the  use  of  contrast  materials,  and  has  been  previously 
demonstrated  in  brain  and  kidney.  The  spin  tagging  technique  may  provide  more 
sensitive  and  specific  assessment  of  tissue  parameters,  at  higher  spatial  resolutions,  by 
virtue  of  its  ability  to  improve  signal  to  noise  ratio  (SNR)  over  a  reasonable  scan  time. 
This  is  another  advantage  over  the  dynamic  contrast  enhancement  technique,  which 
has  a  limited  scan  time  and  thus  a  limited  resolution  and  SNR.  A  limited  scan  time  is 
necessary  to  avoid  the  washout  of  contrast  material. 

Studies  have  shown  that  malignant  tumors  have  higher  water  content  than 
normal  tissues  and  have  higher  perfusion  than  both  normal  tissues  and  benign  lesions. 
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This  is  because  malignant  tumors  induce  high-level  angiogenesis  resulting  in  the 
increase  of  vascularity  [25].  The  rapid  growth  of  malignant  tumors  suggests  high 
metabolic  rates  at  the  tumors,  which  require  an  increased  supply  of  nutrients  and 
oxygen  as  well  as  the  increased  effort  to  remove  the  waste  materials.  Therefore,  a 
higher  blood  perfusion  and  a  higher  blood  volume  at  malignant  tumors  than  normal 
tissue  and  benign  lesions  are  expected  [26,  27,  28].  High  water  content  leads  to  a  long 
MR  spin-lattice  relaxation  time  T\.  During  the  continuous  research  effort  for  this 
dissertation,  the  spin  tagging  technique  that  has  been  developed  [29,30,31,32,33] 
measures  7i  and  relative  blood  perfusion.  These  two  parameters  are  used  as  the 
primary  factors  in  determining  the  malignancy  in  the  breast  cancer  detection  technique 
developed.  This  technique  is  still  in  its  exploration  state,  and  improvement  is  still 

I 

required. 

Microcirculation  in  malignant  tumors  is  complex.  But  macroscopically,  tumor 
vasculature  can  be  studied  in  terms  of  two  idealized  models:  peripheral 
vascularization  where  vessels  are  localized  primarily  at  the  periphery,  and  central 
vascularization  where  vessels  proliferate  from  the  center  like  branches  from  a  tree 
[28].  The  spin  tagging  technique  developed  focused  on  investigating  this  macroscopic 
view  of  the  tumor  microcirculation,  instead  of  the  details  of  the  microcirculation  itself. 

The  rate  of  magnetization  transfer  between  macromolecules  and  water  has  also 
been  shown  to  be  a  valuable  indicator  of  malignant  disorder  in  the  breast  [34,35]  and 
other  parts  of  the  body  [36].  The  level  of  magnetization  transfer  can  be  used  to 
complement  the  application  of  T\  and  perfusion  in  malignancy  detection. 
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Kallinowski's  group  found  that  perfusion  has  a  wide  range  in  malignant  tumors  [26], 
Magnetization  transfer  might  become  more  important  if  perfusion  fails  in  tumor 
recognition. 

Based  on  the  above  three  parameters,  and  additional  statistical  factors,  a  simple 
malignancy  decision  technique  based  on  thresholds  and  a  more  sophisticated  decision 
technique  based  on  Bayesian  theory  are  developed.  The  former  technique  has  been 
implemented  in  a  visualization  program  (called  BreastView)  to  aid  researchers  and 
clinicians  in  decision  making.  The  ROC  (receiver  operating  characteristic)  curve  based 
on  the  thresholds  or  the  risk  ratios  used  in  Bayesian  decision  theory  can  be  used  to 
evaluate  the  technique  presented. 

The  following  six  chapters  (Chapters  4  to  9)  will  be  dedicated  to  the  discussion 
of  the  arterial  spin  tagging  technique  development  and  its  application  to  breast  cancer 
detection.  Chapter  4  describes  the  concept  of  perfusion  and  the  theory  of  the  spin 
tagging  technique.  Chapter  5  discusses  the  arterial  spin  tagging  pulse  sequences 
developed.  They  include  the  fast  SPGR  (SPoiled  Gradient  Recalled)  based  and  EPI 
(Echo  Planer  Imaging)  based  sequences.  Magnetization  transfer  and  its 
implementation  are  also  discussed.  In  Chapter  6,  the  data  processing  techniques, 
including  the  semi-log  linear  regression  technique  and  the  minimum-mean-square- 
error  technique,  are  discussed.  Chapter  7  discusses  the  techniques  of  image 
processing.  Chapter  8  contains  the  results  and  discussion.  Chapter  9  summarizes  this 
project,  presents  conclusions  and  suggests  areas  of  further  development. 
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Chapter  4 


The  Concept  of  Perfusion  and  the  Theory  of  Spin  Tagging 

Before  proceeding  to  the  discussion  of  spin  tagging,  the  concept  of  perfusion 
needs  to  be  clarified.  Buxton  et  al  well  defines  the  concept  of  perfusion  appropriate  for 
imaging  in  Clinical  Magnetic  Resonance  Imaging  edited  by  Edelman  et  al  [37]. 
Considering  a  volume  of  tissue  (V)  supplied  with  arterial  blood  at  a  rate  F  (milliliters 
per  minute)  to  the  capillary  beds  within  V,  the  perfusion  (f)  is  then /=  F/V  (the 
milliliters  of  arterial  blood  delivered  per  minute  per  milliliter  of  tissue).  With  this 
definition  of  perfusion, /has  dimensions  of  1/minute  and  is  thus  a  rate  constant  (Figure 
4.1).  This  definition  differs  slightly  from  the  standard  definition  of  perfusion  as 
milliliters  of  blood  per  minute  per  100  g  of  tissue.  The  latter  definition  is  more 
practical  for  experimental  studies  in  which  tissue  samples  are  measured  for 
concentrations  of  radioactive  tracers,  but  the  former  definition  is  convenient  for 

t 

imaging  studies  where  a  known  volume  of  tissue  is  imaged  [37]. 

In  our  research  protocol,  approximately  10  slices  within  the  breast,  each  with  a 
thickness  of  3  mm,  are  studied.  Each  slice  is  individually  examined  to  investigate 
blood  perfusion.  The  perfusion  is  measured  at  each  voxel  in  these  slices.  In  breast 
tissue,  the  blood  vessels  are  composed  of  capillaries,  connected  by  a  few  larger  vessels 
(arterioles,  metarteriles  and  venules).  The  diameters  of  capillary,  arteriole,  metarterile 
and  venule  are  5  to  10  pm,  5  to  100  pm,  10  to  20  pm  and  5  to  100  pm  respectively. 
Blood  flows  slowly  from  arteriole  to  metarterile  to  the  capillary  bed  and  then  to  the 
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(b) 

Figure  4.1.  The  meaning  of  perfusion.  Two  idealized  capillary  beds  illustrate  the 
distinction  between  blood  volume  (Vb),  blood  velocity  (u)  and  perfusion  (f).  The  Vb’s  as 
well  as  the  u’s  are  the  same  in  Configurations  (a)  and  (b),  but /in  Configuration  (a)  is 
twice  as  large  as  in  Configuration  (b).  In  Configuration  (a),  a  tissue  volume  V  is  fed  by 
two  arteries  each  with  flow  F.  In  Configuration  (b),  the  capillaries  are  twice  as  long  so 
that  the  blood  volume  is  the  same,  but  fed  by  only  one  artery. 
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venules.  Water  molecules  can  move  freely  between  the  capillaries  and  their 
surroundings.  The  average  velocity  of  blood  flow  in  capillaries  is  1  mm/sec  [2],  The 
voxel  resolution  used  in  our  research  protocol  is  1 .33  mm  x  1 .33  mm  x  3  mm,  which 
is  also  the  range  used  in  typical  clinical  MR  imaging.  There  are  about  50  capillaries 
penetrating  each  voxel  in  various  directions.  Therefore,  perfusion  represents  the 
combined  activity  of  micro  flows  within  a  voxel  (Figure  4.2).  A  two-component 
model  is  used  to  approximate  the  activity  of  the  water  molecules.  In  this  model,  water 
molecules  in  blood  vessels  entering  a  slice  of  breast  tissue  mix  thoroughly  with  water 
molecules  outside  the  vessels  before  they  leave  the  slice.  For  this  two-component 
model,  the  Bloch  Equation  for  the  breast  tissue  water  magnetization  at  the  slice  of 
tissue  investigated  after  any  RF  perturbation  can  be  written  as  [24] 


dM(t)  M0  -  M(t) 
dt  Tx 


-  kforM(t)  +  krev  Mm(t)  +fMa(t)  -  f 


M(t) 

A  ) 


(4.1) 


where, 

* 

M(t )  =  the  longitudinal  water  magnetization  of  a  voxel  of  breast  tissue  at  time 


t, 

Mo  =  the  value  of  M(t )  under  fully  relaxed  conditions, 

T\  =  the  spin-lattice  relaxation  time  constant  of  the  breast  tissue, 

Mm(t)  =  the  longitudinal  magnetization  of  macromolecules  in  a  voxel  of  tissue, 
kfor ,  krev  =  the  magnetization  transfer  rate  constants  from  water  to 
macromolecules  and  vice  versa, 
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The  slice  of  tissue 
to  be  investigated 
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be  investigated.  A  voxel  within  the  slice  is  identified.  As  shown,  small  blood  vessels  randomly 
penetrate  the  voxel  resulting  in  random  directions  of  blood  flow  (indicated  by  long  arrows).  Water 
molecules  (indicated  by  short  arrows)  move  easily  in  and  out  of  the  blood  vessels.  The  blood 
perfusion  of  the  voxel  represents  the  overall  activity  of  the  micro  flows  within  it. 


Ma  =  the  longitudinal  water  magnetization  per  ml  of  arterial  blood  entering  a 
voxel  from  outside  of  the  slice  investigated, 

/=  perfusion  in  ml  of  arterial  blood  per  second  per  voxel  (with  dimensions  of 
sec'1), 

X  =  the  ratio  of  water  content  between  breast  tissue  overall  and  the  arterial 
blood  within  it. 

To  simplify  Eq.  (4.1),  we  need  to  make  two  assumptions.  The  first  assumption 
is  that  the  effect  of  magnetization  transfer  is  negligible,  specifically. 

Then  Eq.  (4.1)  can  be  simplified  to 

=  M,  -M(f)  +  (,)  .  j  (4.2) 

Eliminating  the  magnetization  transfer  terms  can  also  be  rationalized  as  being 
included  in  T\  relaxation.  The  effect  of  magnetization  transfer  has  to  be  measured 
using  other  means  to  improve  the  precision  of  Eq.  (4.1),  such  as  the  use  of  binomial 
pulses  [38]  or  off-center  RF  irradiation  [36,39]. 

Various  spin  tagging  techniques  have  been  used  to  estimate  perfusion /  (or 
/ 12).  Some  are  qualitative  [40]  and  others  are  more  quantitative  [24,41,42].  The 
various  techniques  are  all  based  on  three  different  spin  states:  equilibrium,  saturation 
and  inversion. 
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The  arterial  spin  tagging  method  presented  in  this  dissertation  is  also  based  on 
two  of  these  three  states:  the  equilibrium  and  inversion  states,  which  should  give  the 
maximum  spin  signal  contrast.  The  measurement  of  both  T\  and / // 1  are  achieved.  This 
general  scheme  of  the  pulse  sequence  is  based  on  inversion  recovery.  Spins  in  a  slab 
of  breast  tissue  are  inverted  first  using  a  RF  pulse.  The  thickness  of  the  slab  can  be 
controlled  by  the  strength  of  the  slab  selective  gradient  (normally  a  Z  gradient).  Then 
the  levels  of  recovery  to  equilibrium  are  measured  at  different  times  (Figure  4.3).  The 
rate  of  recovery  can  be  estimated  from  the  recovery  curve. 

Parameters  T\  and / IX  in  Eq.  (4.2)  can  be  estimated  from  two  spin  tagging 
conditions.  Under  Condition  1  (Selective  tagging  condition),  spins  in  the  slice  of 
breast  tissue  to  be  investigated  are  inverted  first,  and  then  the  signal  corresponding  to 

1 

the  overall  spin  magnetization  from  the  same  slice  is  acquired  at  a  later  time,  which  is 
varied  to  obtain  a  good  spin  recovery  curve.  Since  the  spins  entering  the  slice  of  breast 

tissue  are  not  perturbed,  they  are  in  equilibrium.  Then, 

* 

Ma(t)  =  Ma 

where, 

Ma  =  the  value  of  Ma(t)  under  the  fully  relaxed  condition. 

In  the  two-component  model,  we  assume  that  the  spins  in  arterial  blood 
entering  the  slice  can  mix  thoroughly  with  spins  in  the  surrounding  breast  tissue  before 
they  leaving  the  slice.  Thus  in  the  equilibrium  condition, 
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Signal  recovery  after  time  77 


77 


Figure  4.3.  Signal  recovery  after  spin  inversion.  After  the  spin  dipoles  are 
inverted,  they  relax  back  to  their  equilibrium  positions  in  an  exponential 
fashion.  Shown  here  is  the  signal  acquired  along  a  typical  spin  recovery 
curve.  An  image  only  shows  the  magnitude  of  the  actual  signal  measured. 
The  rate  of  recovery  can  be  determined  by  the  signal  intensity  in  images 
acquired  at  different  points  of  the  recovery  curve.  77  is  a  point  in  time  at 
which  an  image  is  acquired. 


Eq.  (4.2)  then  becomes 


fMl  =  f 


Eq.  (4.3)  is  the  standard  relaxation  equation  below 

dM(t)  M0  -  M(t) 


dt 


lu 


with 


(4.3) 


Tu.  '  Ti  * 


(4.4) 


where  Tis  is  the  apparent  spin-lattice  relaxation  time  measured  under  the  Selective 
condition.  The  subscript  "s"  means  "under  Selective  condition". 

Under  Condition  2  (Non-selective  tagging  condition),  spins  are  inverted 
everywhere,  but  signals  are  acquired  only  from  the  slice  to  be  investigated.  The  spins 
inside  and  outside  the  slice  of  breast  tissue  are  thus  the  same.  Then 


and  Eq.  (4.2)  becomes 


fMa{t)  =  / 


Mjt) 

A, 


dM(t )  M0  -  M(t) 

dt  t; 
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The  T\  measured  in  this  case  is  the  true  T\  without  the  influence  of  blood  flow,  which 


can  be  denoted  as 


Tm  =  Tx  (4.5) 

and 


The  subscript  "n"  denotes  the  T\  and  the  b  value  measured  under  the  Non-selective 
condition,  distinguishing  them  from  the  Selective  case. 

Based  on  the  bs  and  b„  values  obtained  from  the  two  tagging  conditions,  we 
can  determine  the  perfusion  parameter 

J  =  bs  -  K  (4.7) 

V 

The  term / /A  can  be  interpreted  as  the  rate  of  spin  replacement  in  a  voxel  in  the 
direction  perpendipular  to  the  slice.  Compared  to  normal  breast  tissue  and  benign 
lesions,  we  would  expect  a  higher  / 1  A  in  malignant  tumors. 

In  summary,  this  chapter  has  covered  the  concept  of  perfusion  and  its 
measurement  in  MRI.  The  theory  of  spin  tagging  has  been  thoroughly  discussed,  as 
well  as  how  it  is  used  to  measure  the  perfusion  parameter / 7 A.  The  next  chapter  will 
cover  the  arterial  spin  tagging  pulse  sequences  developed  for  the  measurement  of  the 
TVs  under  the  Selective  and  Non-selective  conditions,  from  which  the  perfusion 
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parameter/ IX  can  be  calculated,  and  for  the  measurement  of  the  level  of 
magnetization  transfer  between  macromolecules  and  water. 
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Chapter  5 


Arterial  Spin  Tagging  Pulse  Sequences 

To  capture  the  ideal  spin  recovery  after  the  spin  tagging  inversion  pulse,  some 
design  criteria  need  to  be  considered:  (1)  The  slice  profiles  generated  by  the  inversion 
RF  pulse  and  the  data  acquisition  RF  pulse  should  match  perfectly  with  each  other 
during  the  Selective  condition  so  that  the  inverted  spins  are  the  spins  being  measured. 
(2)  The  total  time  of  whole-image  data  acquisition  for  each  point  of  the  recovery  curve 
should  be  as  short  as  possible  to  avoid  the  averaging  effect.  (3)  The  total  scan  time  for 
the  full  breast  examination  should  be  reasonably  short.  (4)  The  data  acquisition  should 
not  disturb  the  natural  spin  recovery  after  the  inversion  RF  pulse  excitation. 

To  satisfy  Criterion  1,  an  optimized  RF  inversion  pulse  and  an  optimized  RF 
acquisition  pulse  have  been  designed  based  on  the  Shinnar  Le-Roux  (SLR)  algorithm 

[43,12,15].  These  optimized  RF  pulses  allow  sharp  rectangular  slice  profiles  so  that 

* 

the  spins  being  tagged  are  basically  the  spins  being  measured  under  the  Selective 
condition  (Figure  5.1).  The  section  of  inverted  magnetization  (slab)  is  created  by 
applying  the  SLR  inversion  pulse  in  the  presence  of  the  slice  selection  gradient.  The 
thickness  of  the  inversion  slab  can  be  controlled  by  changing  the  gradient  strength.  By 
applying  the  SLR  pulse  without  an  associated  slice  selection  gradient,  a  non-selective 
slab  containing  inverted  spins  in  the  entire  sensitive  volume  of  the  transmit  RF  coil  is 
created.  To  optimize  the  RF  pulse  profiles,  tradeoffs  have  to  be  considered  among  the 
length,  the  amplitude  and  the  bandwidth  of  the  RF  pulses.  In  general,  a  better  slice 
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Inversion  RF  Amplitude 


Time  (ms) 

(a) 

10°  Acquisition  RF  amplitude 


Time  (ms) 

(c) 


Mz  Amplitude 


Frequency  (kHz) 

(b) 

Mz  Magnitude  &  Phase  response 


Figure  5.1.  SLR  inversion  and  data  acquisition  RF  pulses,  (a)  and  (b)  show  the 
SLR  inversion  RF  pulse  and  the  near  rectangular  frequency  response  for  the 
longitudinal  magnetization,  (c)  shows  the  10°  SLR  data  acquisition  RF  pulse, 
(d)  shows  the  magnitude  of  the  frequency  response  of  the  longitudinal 
magnetization  and  the  associated  phase. 
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profile  is  obtained  by  using  a  RF  pulse  with  a  longer  pulse  width,  a  higher  amplitude 
and  a  larger  bandwidth.  However,  a  long  data  acquisition  RF  pulse  can  lead  to  a 
lengthening  of  the  scan  time  because  it  is  used  many  times  for  the  completion  of  the 
data  acquisition  for  the  whole  image.  Its  excitation  effect  on  the  moving  spins  due  to 
perfusion  might  also  change  during  a  long  RF  excitation.  The  high  amplitude  of  a  RF 
pulse  can  go  beyond  the  limitation  of  the  hardware.  For  a  fixed  pulse  length,  a  larger 
bandwidth  can  lead  to  a  better  slice  profile,  but  it  also  induces  a  higher  RF  amplitude. 
There  is  a  hardware  limitation  on  the  bandwidth  itself,  but  it  has  not  led  to  any 
problem  for  the  RF  pulse  design  in  this  dissertation.  The  final  design  of  the  RF  pulses 
is  shown  on  Table  5.1  below: 


Table  ‘5.1 .  The  SLR  RF  Pulse  Design  Specifications 


Parameters 

Inversion  Pulse 

Data  Acquisition  Pulse 

Pulse  Type 

Inversion 

Spin  echo 

Phase  Type 

Maximum  phase 

Linear  phase 

Pulse  duration  (ms) 

15.6 

6.4 

Pulse  bandwidth  (kHz) 

0.95 

1.25 

Ripple  in  passband 

1.50% 

5% 

Ripple  in  stopband 

0.40% 

5% 

Resolution  (#  of  points) 

1560 

640 

Flip  angle  (degree) 

180 

10 

The  above  design  specifications  were  implemented  using  the  RF  design  tool  that 
comes  with  the  EPIC  software.  The  resulting  attributes  for  the  RF  pulses  are  listed  in 
Table  5.2  below: 
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Table  5.2.  The  RF  Resulting  Attributes  Based  on 
the  Design  Specifications  in  Table  5.1 


Parameters 

Inversion  Pulse 

Data  Acquisition  Pulse 

Actual  Resolution  (#  of  points) 

1558 

638 

isoDelay  (ms) 

unknown 

3.195 

Pulse  area  (G  ms) 

0.0343 

0.1216 

Absolute  width  (ratio) 

0.1091 

0.2309 

Effective  width  (ratio) 

0.0445 

0.1249 

Duty  cycle 

13.22% 

24.14% 

Maximum  pulse  width 

8.92% 

24.14% 

Maximum  B1  (G) 

0.219366 

0.00838565 

Maximum  integral  Bl2  (G2) 

0.0334052 

5.62188  x  10‘5 

Maximum  RMS  Bl  (G) 

0.0462748 

0.00296381 

IsoDelay:  Time  in  ps  from  the  zero  phase  reference  of  the  pulse  to  the  end  of 


pulse. 


Pulse  area:  The  integrated  area  of  the  pulse. 

Absolute  width:  The  ratio  of  absolute  pulse  area,  calculated  by  integrating  the 

absolute  value  of  pulse  amplitude  pver  time,  to  that  of  a  hard  pulse  of  unit  amplitude 

i 

and  equal  duration. 

Hard  pulse:  A  rectangular  RF  pulse,  also  called  a  standard  RF  pulse. 

Effective  width:  Ratio  of  pulse  energy  to  that  of  a  hard  pulse  of  unit  amplitude 
and  equal  duration. 

Duty  cycle:  Duty  cycle  of  pulse  defined  by  GE  Medical  Systems  Inc..  It  is  the 
percentage  of  pulse  duration  in  which  amplitude  exceeds  22.3%  of  maximum  pulse 
amplitude. 
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Maximum  pulse  width:  Percentage  of  pulse  width  that  exceeds  22.4%  of 
pulse's  maximum  amplitude  that  is  found  in  the  largest  lobe  of  the  pulse. 

B1 :  The  magnetic  field  of  the  RF  pulse. 

Maximum  B1 :  Peak  B1  field  strength  in  Gauss. 

Maximum  integral  B1  :  Peak  integral  of  B1  in  Gauss  . 

Maximum  RMS  Bl:  Peak  RMS  B1  in  Gauss. 

The  above  attributes  are  used  to  scale  the  amplitudes  of  the  RF  pulses  against 
the  actual  amplitude  of  a  90°  sine  RF  pulse  that  causes  a  true  90°  spin  flip  determined 
in  a  prescan.  They  are  also  used  to  calculate  the  heat  absorption  rate  to  the  patient  so 
that  a  maximum  repetition  rate  of  the  RF  pulse  excitation  can  be  set  according  to  the 
safety  guidelines  that  have  been  established  by  health  agencies  such  as  the  FDA. 

i 

The  RF  pulses  designed  might  not  be  optimized  in  terms  of  the  flip  angle  and 
pulse  width.  They  can  be  manipulated  without  changing  the  shape  of  the  RF  pulse  or 
the  slice  profile.  But  the  peak  amplitude  of  the  RF  pulse  is  changed  according  to 


Bl 


peak 


=  Bl 


nominal 

peak 


flip  nom_  pw 

nom_  flip  pw 


(5.1) 


where, 


flip  =  the  flip  angle  of  the  RF  pulse  used, 

pw  =  the  pulse  width  of  the  RF  pulse  used, 

nom  Jlip  =  the  normal  flip  angle  of  the  RF  pulse  designed, 

nom _pw  =  the  normal  pulse  width  of  the  RF  pulse  designed, 

Blpeak  =  the  peak  amplitude  of  the  RF  pulse  used, 
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B l^mal  =  the  peak  amplitude  of  RF  pulse  at  normal  flip  angle  and  pulse  width. 

Normally,  only  the  flip  angle  is  manipulated  during  scanning  because  the  change  of 
the  pulse  width  leads  to  a  change  of  the  timing  diagram.  For  the  spin  tagging  sequence 
implemented  here,  the  inversion  pulse  is  fixed.  The  optimal  flip  angle  for  the  data 
acquisition  RF  pulse  has  been  determined  by  inspecting  the  resulting  images  at  various 
flip  angles. 

To  satisfy  Criteria  2  and  3,  fast  data  acquisition  techniques  should  be  used. 
These  techniques  can  include  fast  SPGR  (SPoiled  Gradient  Recalled),  EPI  (Echo 
Planar  Image)  and  Spiral  sequences.  They  all  have  their  own  technical  advantages  and 
disadvantages.  The  fast  SPGR  sequence  has  been  tested  repeatedly  in  both  research 

I 

and  clinical  environments.  The  image  quality  is  good  and  the  images  contain  few 

artifacts.  It  is  much  less  demanding  on  the  gradient  amplitude  and  the  receiver 

bandwidth  than  the  other  two  sequences.  However,  it  is  much  slower  than  the  other 
* 

two  techniques  and  thus  is  more  prone  to  motion  artifacts.  It  also  has  more  chance  of 
disturbing  the  natural  spin  recovery  due  to  the  need  of  repeated  application  of  the  data 
acquisition  RF  pulse  for  the  whole  image.  On  the  other  hand,  the  EPI  and  Spiral 
sequences  only  require  one  application  of  the  data  acquisition  RF  pulse  for  the  whole 
image  or  for  many  ky  lines  in  the  case  of  interleaved  data  acquisition.  They  are  more 
capable  of  preserving  the  natural  spin  recovery  curve.  They  also  reduce  scan  time  and 
thus  are  less  prone  to  motion  artifacts.  However,  they  demand  a  high  amplitude  and  a 
fast  slew  rate  that  some  MR  systems  cannot  provide.  Due  to  the  use  of  fast  alternating 
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gradient  pulses,  they  also  lead  to  serious  artifacts  and  a  lower  image  quality.  EPI 
appears  advantageous  compared  to  Spiral  imaging.  An  EPI  sequence  is  easier  to 
implement,  allows  for  more  control  over  the  frequency  and  phase  encoding  directions, 
and  is  capable  of  asymmetric  matrix  data  acquisition.  The  techniques  for  the  EPI 
artifact  removal  and  interleaved  data  acquisition  have  been  folly  developed.  Although 
the  Spiral  sequence  is  more  efficient  in  the  use  of  the  gradients,  it  is  also  more 
demanding  on  the  gradients. 

The  Signa  Advantage  1.5  T  GE  MR  system  (GE  Medical  Systems  Inc., 
Milwaukee)  has  been  the  system  available  in  the  research  site  at  University  of 
California  Davis  Medical  Center  during  the  development  of  this  project.  This  system 
could  not  support  the  EPI  and  spiral  application  for  this  project,  so  the  arterial  spin 
tagging  sequence  has  been  developed  based  on  a  fast  SPGR  data  acquisition  scheme. 
This  spin  tagging  sequence  has  been  fullyf  ested  on  phantoms,  normal  subjects  and 
patients,  and  has  been  proven  to  function  as  expected.  The  data  used  for  analysis  is 
mainly  from  the  application  of  this  sequence.  In  April  1999,  the  Signa  Advantage 
system  was  upgraded  to  Signa  Horizon  LX  1.5  T  GE  MR  system  (GE  Medical 
Systems  Inc.,  Milwaukee).  The  spin  tagging  sequence  based  on  the  fast  SPGR  data 
acquisition  scheme  has  been  converted  to  run  on  this  new  system.  But  this  new  system 
still  has  problems  supporting  the  polarity  alternation  of  the  EPI  X  gradient  pulse  train 
from  one  interleave  to  the  next.  Therefore,  the  EPI  based  arterial  spin  tagging  pulse 
sequence  has  only  been  developed  to  run  on  computer  simulation. 
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The  180°  RF  refocusing  pulse  has  been  eliminated  from  all  of  the  three  fast 
sequences  to  achieve  fast  data  acquisition.  Since  each  data  acquisition  period  is  very 
short,  the  transverse  magnetization  is  not  completely  dephased.  This  problem  is  most 
severe  for  fast  SPGR  because  its  TR  period  is  15  ms  as  compared  to  about  100  ms  for 
EPI  or  Spiral.  To  avoid  the  transverse  magnetization  being  carried  over  to  the  next 
data  acquisition  periods,  it  needs  to  be  dephased  or  randomized.  This  has  been 
achieved  by  using  a  combination  of  gradient  and  RF  spoiling  [44].  The  spoiling 
techniques  will  be  investigated  by  looking  into  the  fast  SPGR  pulse  sequence 
following  (Figure  5.2). 

5.1.  Arterial  Spin  Tagging  Sequence  Based  on  a  Fast  SPGR  Data  Acquisition 
Scheme 

1 

The  fast  SPGR  sequence  is  a  gradient  recalled  echo  (GRE)  based  sequence,  so 
it  looks  similar  to  the  one  discussed  in  Chapter  2  (Figure  2.6).  The  optional  killer  pulse 
in  the  X  gradient  has  not  been  used.  To  achieve  a  short  repetition  period,  a  higher  data 
acquisition  frequency  is  normally  used.  This  leads  to  a  higher  read-out  gradient 
amplitude  (the  gxw  pulse  amplitude).  Shorter  pulse  widths  are  also  used  in  other 
gradients,  except  the  slice-selection  gradient  that  is  determined  by  the  slice  thickness 
and  the  bandwidth  of  the  RF  pulse.  The  other  method  used  to  shorten  the  TR  period  is 
the  concurrent  application  of  the  gradient  pulses  except  during  RF  pulse  excitation  and 
data  acquisition.  This  fast  SPGR  sequence  has  been  implemented  for  clinical  studies 
by  GE  Medical  Systems  Inc.,  and  the  image  quality  from  this  sequence  has  also  been 
verified  [44]. 
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TR  (time  of  repetition) 


Figure  5.2.  A  fast  SPGR  sequence.  The  gradient  and  RF  pulses  played  in 
one  repetition  period  of  a  fast  SPoiled  Gradient  Recalled  (SPGR)  sequence. 
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The  details  for  spoiling  gradient  pulse  design  are  described  by  Zur  and  his 
colleagues.  They  suggested  that  a  phase  increment  of  1 17°  or  123°  at  each  TR  period 
would  give  adequate  cancellation  of  transverse  magnetization  in  clinical  MRI  [9].  In 
the  SPGR  sequence  implemented  by  GE  Medical  Systems  Inc.,  the  pulse  area  of  the 
spoiling  gradient  pulse  has  been  designed  to  have  at  least  4 n  of  phase  dispersion 
across  the  slice  at  each  TR  period,  combined  with  the  net  phase  introduced  in  the 
slice-selection  gradient  pulse  and  its  rewinding  pulse.  But  the  pulse  area  range  is  set 
between  100  ms •  G/cm  and  400  ms*  G/cm.  The  pulse  area  calculation  is  given  by  the 
following  equation: 


area_  gzk  =  — - (/_  exa  -  ?w-^zr^—) .  |  a_  gzrf\\ 

y-Az  2 


where, 


(5.2) 


area_gzk  =  the  pulse  area  of  the  spoiling  gradient  pulse  (sec  •  G/cm), 

Az  =  the  slice  thickness  (cm), 

* 

y-  the  gyromagnetic  ratio  for  the  proton  =  4258  Hz/G, 
a_gzrfl  =  the  amplitude  of  the  gzrfl  pulse  (G/cm), 
pw_gzrfla  =  the  width  of  the  ramp  time  of  the  gzrfl  trapezoid  pulse  (sec), 
t_exa  =  the  time  from  the  start  to  the  middle  of  RF  pulse,  including  pw_gzrfla 
(sec). 

To  achieve  effective  spoiling  after  just  a  few  or  no  TR  periods,  RF  spoiling  is 
also  implemented.  For  the  SPGR  sequence  implemented  by  GE  Medical  Systems  Inc., 
the  RF  spoiling  use  the  following  scheme: 
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(1)  Initialization: 

fs_pi  =  phase  board  value  for  n  radians  =  32752, 
fs_2pi  =  phase  board  value  for  2n  radians  =  65504, 
seq  =  the  fast  SPGR  repetition  number  =  0, 
seed  =  21001. 

(2)  Looping  from  one  TR  period  to  the  next: 

phase  =  (phase  +  seq  *  seed  +  3  *  fs_pi)  %  fs_2pi  -  fs_pi; 
seq  ++; 

An  application  of  an  inversion  RF  pulse  before  the  full  data  acquisition  scheme 
is  the  basis  of  the  arterial  spin  tagging  sequence.  An  inversion  RF  pulse  is  applied  first 
to  tag  the  spins.  Data  is  then  acquired  for  the  entire  image  using  a  fast  SPGR  data 
acquisition  scheme.  In  the  standard  clinical  sequence,  data  acquisition  starts  at  one  end 
of  the  A  space  and  a  significant  amount  of  time  elapses  before  Ay  =  0  occurs.  This  time 
is  too  long  to  map  out  the  T\  recovery  curve.  The  original  clinical  sequence  must  be 
modified  to  support  a  symmetric-A  space  data  acquisition  (Figure  5.3).  Data 
acquisition  starts  at  the  Ay  =  0  line  of  the  A  space,  and  then  alternates  below  and  above 
this  Ay  line  until  the  lowest  and  highest  Ay  lines.  This  allows  the  signal  of  the 
reconstructed  imdge  to  have  the  highest  weighting  on  the  data  points  at  the  beginning 
of  the  acquisition.  Thus,  the  effective  time  of  inversion,  which  is  the  time  from  spin 
inversion  to  the  effective  zero  A  space,  can  have  a  wide  range  and  thus  map  the  T\ 
recovery  curve.  Since  the  RF  phase  is  alternated  in  fast  SPGR,  the  order  of  data 
acquisition  has  been  constructed  in  a  way  such  that  this  alternating  RF  phase  pattern  is 
retained  in  both  the  actual  data  acquisition  as  well  as  the  whole  data  set.  This  allows 
for  the  correct  image  reconstruction  using  the  algorithms  installed  in  the  MR  system. 
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Figure  5.3.  The  k  space  for  the  fast  SPGR  symmetric  data 
acquisition.  The  total  number  of  ky  lines  is  Nmax. 


The  complete  fast  SPGR  based  arterial  spin  tagging  pulse  sequence  is  shown  in  Figure 
5.4. 

Odd-number  interleaved  &-space  acquisition  [45]  has  also  been  implemented  to 
allow  high  spatial  resolution  in  the  phase-encoding  direction.  The  k  space  is  completed 
from  different  inversion  periods  (interleaves).  By  using  an  odd  number  of  interleaves, 
the  data  combined  would  appear  obtained  from  a  single-shot  acquisition.  This  would 
allow  the  simple  image  reconstruction  as  before.  But  the  scan  time  would  increase 
with  the  number  of  interleaves.  Since  then,  a  more  effective  technique  to  approach  a 
sufficiently  high  resolution  has  been  found  by  using  the  half  phase  field  of  view  data 
acquisition.  In  this  approach,  the  field  of  view  (FOV)  in  the  Y  direction  is  only  half  of 
that  in  the  X  direction.  This  takes  advantage  of  the  fact  that  two  breasts  are  scanned  at 

i 

a  time.  The  long  axis  can  be  chosen  as  the  frequency  direction  (X  direction).  The  FOV 
required  in  the  Y  direction  is  obviously  only  half  of  that  in  the  X  direction.  So  only 
half  of  the  phase-encoding  steps  are  needed  to  achieve  the  same  spatial  resolution  as 
before.  This  allows  the  doubling  of  the  spatial  resolution  without  any  sacrifice  of  the 
scan  time.  Certainly,  if  the  resolution  is  not  sufficient  after  the  application  of  this 
technique,  the  odd-number  interleaved  A>space  acquisition  can  be  added  on  to  it.  Or  if 
the  geometry  of  the  object  does  not  allow  the  half-phase  FOV  data  acquisition 
technique,  the  odd-number  interleaved  /c-space  acquisition  technique  can  be  used. 

For  this  fast  SPGR  spin  tagging'  sequence  operated  in  the  Signa  Advantage  MR 
system,  the  following  control  variables  have  been  used  for  all  our  studies: 

TE  =  9.244  ms  (The  minimum). 
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optr  =  the  time  of  repetition  for  individual  data  acquisition  frame  =  16.62  ms 
(The  minimum), 

Receiving  bandwidth  =  31.25  kHz, 

Data  acquisition  flip  angle  =  10°, 

Slice  thickness  =  3  mm, 

Field  of  view  (FOV)  in  X  direction  =  34  cm, 

Field  of  view  (FOV)  in  Y  direction  =  17  cm  (by  using  half-phase  FOV), 

X  matrix  resolution  =  256, 

Y  matrix  resolution  =  120, 

Number  of  interleaves  =  1  (single-shot), 

Time  of  inversion  =  start  with  615  ms,  decrement  to  15  ms  with  a  step  size  of 

i 

100  ms. 

The  total  scan  time  for  the  acquisition  of  data  for  the  entire  image  (TR  in 

Figure  5.4)  is  2.7  sec,  including  the  inversion  RF  pulse.  The  data  acquisition  under  the 

i 

Regular  condition  is  designed  simply  by  setting  the  amplitude  of  the  inversion  RF 
pulse  and  the  associated  Z  gradient  pulse  to  zero,  and  thus  keeping  the  same  timing  for 
the  gradient  and  RF  pulses. 

The  Signa  Horizon  LX  MR  system  has  a  higher  gradient  amplitude  tolerance. 

A  fast  data  acquisition,  which  requires  a  higher  read-out  gradient  amplitude,  is  used, 
and  thus  the  data  acquisition  period  (optr)  can  be  shortened  to  12.5  ms.  The  TE  now 
becomes  5.2  ms.  Since  the  TR  is  shortened,  a  higher  Y  matrix  resolution  can  be  used. 
The  default  Y  matrix  resolution  is  set  at  128.  The  total  time  for  the  whole  image 
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acquisition,  including  the  inversion  RF  pulse,  is  dropped  down  to  2.22  sec.  This  TR 
period  should  be  designed  within  a  reasonable  range  (2  sec  to  3  sec  for  tissue)  so  that 
sufficient  signal  can  be  recovered  to  have  a  good  SNR.  The  rest  of  the  control 
variables  remain  the  same. 

For  each  breast  study,  approximately  10  equally  spaced  slices  in  the  breasts  are 
imaged  with  a  total  scan  time  of  approximately  25  minutes.  Two  breasts  are  scanned 
simultaneously.  At  each  slice  location,  61  images  are  acquired  in  a  continuous  manner 
with  a  repetition  time  of  2.7  sec  (or  2.22  sec  in  the  LX  system)  based  on  the  arterial 
spin  tagging  sequence  (Figure  5.5).  The  image  acquisition  is  based  on  the  three 
conditions:  Selective,  Non-selective  and  Regular  conditions.  Under  either  the 
Selective  or  Non-selective  condition,  spins  are  inverted  first  and  then  data  acquisition 

i 

starts  at  a  certain  points  (TIbr  s)  of  the  recovery  curve.  For  the  Selective  condition,  only 

spins  within  the  slice  being  measured  are  inverted.  For  the  Non-selective  condition, 

spins  within  the  entire  sensitive  volume  of  the  transmit  RF  coil  are  inverted.  For  either 
* 

Selective  or  Non-selective  condition,  seven  different  time  delays  (or  seven  different 
Tlbr's)  are  used.  Under  each  timing  condition,  four  repeating  images  under  the  same 
constraints  are  acquired.  The  last  three  images  of  each  timing  condition  are  nearly  at 
steady  state  and  are  extracted  for  analysis.  Under  the  Regular  condition,  five  repeated 
images  are  acquired  under  the  same  constraints.  The  last  three  images  are  nearly  at 
steady  state  and  their  mean  is  used  for  analysis.  Using  the  mean  of  three  images 
instead  of  one  helps  to  improve  the  signal-to-noise  ratio  (SNR).  Through  image 
comparison  between  the  Regular  condition  and  the  Selective  or  Non-selective 
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Image  series  plot 


Figure  5.5.  Image  series  plot  using  the  fast  SPGR  based  breast  sequence.  This 
plot  shows  the  signal  at  the  same  pixel  location  across  the  61  images  acquired 
using  this  sequence.  There  are  three  imaging  conditions:  Selective,  Non-selective 
and  Regular  conditions.  For  either  Selective  or  Non-selective  condition,  seven  77 
timing  conditions  are  used.  Four  repeated  images  are  acquired  at  each  77  timing 
condition.  Images  are  acquired  in  a  continuous  manner  with  a  repetition  time  of 
2.7  seconds. 
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condition,  T\s  or  Tin  can  be  calculated  respectively.  The  technique  of  calculation  and 
derivation  is  discussed  fully  in  Chapter  6. 

Appendix  A  outlines  the  detail  of  the  simulation  and  scanner  operation  for  this 
SPGR  based  arterial  spin  tagging  pulse  sequence. 

5.2.  Arterial  Spin  Tagging  Sequences  Based  on  EPI  Data  Acquisition  Scheme 
In  EPI,  the  full-image  data  acquisition  can  be  completed  after  one  RF 
excitation  (Figure  5.6).  EPI  is  still  a  derivative  of  a  GRE  sequence  (Figure  2.6), 
although  they  appear  quite  different.  A  net  zero  phase  accumulation  in  the  X  direction 
is  still  at  the  center  of  each  X  read-out  pulse.  A  set-up  pulse  with  an  area  equal  to  half 
of  the  read-out  pulse  is  applied  before  the  train  of  read-out  pulses.  A  Y-gradient  set-up 
pulse  is  applied  before  the  train  of  X  read-out  pulses.  This  allows  data  acquisition  to 
start  at  the  top  of  the  ky  line.  Y-gradient  blip  pulses  with  an  equal  area  are  applied 
between  the  X  read-out  pulses.  This  allows  the  uniform  decrementation  of  the  ky  lines 

from  the  top  of  k  space  to  the  bottom  of  k  space.  The  net  zero  phase  accumulation  in 

i 

the  Y  direction  occurs  after  half  the  data  has  been  acquired.  Thus,  the  net  zero  phase 
accumulation  considering  both  the  X  and  Y  directions  occurs  at  the  center  of  the  EPI 
train.  Figure  5.7  shows  the  k  space  of  an  EPI  sequence  with  a  matrix  resolution  of  64  x 
64. 

However,  there  is  a  limitation  on  image  matrix  resolution  using  this  one-shot 
technique  because  of  the  T2  decay  and  the  hardware  limitations  on  the  number  of  data 
points  that  can  be  acquired  continuously.  To  achieve  a  good  signal  to  noise  ratio,  all 
the  data  must  be  acquired  before  too  much  signal  is  lost  due  to  T2  decay  (or  T2*  decay 
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Figure  5.6.  (continued) 


* 


Figure  5.7.  The  k  space  for  a  regular  64  x  64  EPI.  This  drawing 
shows  the  k  space  of  the  regular  64  x  64  EPI.  The  direction  of  data 
acquisition  is  from  the  top  to  the  bottom,  and  from  the  left  to  the 
right  and  then  from  the  right  to  the  left  and  so  on. 
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if  field  inhomogeneity  exists).  The  number  of  data  points  that  can  be  acquired 
continuously  is  limited  by  the  receiver.  The  N/2  ghost  also  becomes  more  severe  with 
the  lengthening  of  the  EPI  gradient  pulse  train.  In  his  master's  thesis  the  author  has 
developed  an  odd-number  hybrid  (or  interleave)  EPI  technique  for  high-resolution 
imaging  [45,46,47,48,49].  The  data  for  a  full  k  space  is  completed  using  an  odd- 
number  of  interleaves  (Figure  5.8).  A  time-shift  is  introduced  incrementally  to  the 
interleaves  (Figure  5.9).  This  interleave  technique  allows  the  combined  data  set  to 
appear  that  it  has  been  acquired  in  a  one-shot  EPI.  The  N/2  ghost  appears  the  same  as 
a  one-shot  EPI,  and  thus  can  be  removed  using  the  established  N/2  ghost-removal 
algorithm  [50]. 

To  implement  the  spin  tagging  with  interleaved  EPI,  the  following  steps  are 
performed  at  each  interleave:  (1)  inversion  pulse  and  (2)  20  EPI  acquisitions,  each  of 
which  starts  with  a  10°  RF  pulse  (Figure  5.10).  Each  of  these  20  acquisitions,  with  a 
repetition  time  ( optr  in  Figure  5.10)  of  100  ms,  defines  one  point  in  the  T\  recovery 
curve.  The  total  time  for  each  interleave  (77?  in  Figure  5.10)  is  2  seconds.  The  data  in 
the  other  interleaves  is  gathered  in  the  same  fashion.  Since  repeated  RF  pulses  are 
applied  to  the  same  spins  at  the  same  slice  location,  transverse  magnetization  might 
not  have  completely  dephased  yet.  RF  spoiling  should  be  added  to  the  EPI  sequence 
also.  Figure  5.1 1  shows  the  looping  technique  used  in  the  sequence. 

5.3.  The  Implementation  of  Magnetization  Transfer 

Magnetization  transfer  (MT)  is  only  implemented  for  the  EPI  based  spin 
tagging  sequence.  It  is  not  applicable  to  the  fast  SPGR  based  spin  tagging  sequence. 


Pg.  68 


Figure  5.8.  The  k  space  for  a  128  x  128  odd-number  hybrid  EPI.  This 
drawing  shows  the  k  space  of  the  128  x  128  odd-number  hybrid  EPI 
that  was  developed.  Three  128  x  44  hybrid  scans  are  used.  The  starting 
direction  of  data  acquisition  for  one  hybrid  scan  is  the  opposite  of  the 
previous  one.  The  data  acquisition  direction  of  odd-number  hybrid 
scans  is  from  the  top  to  the  bottom,  and  from  the  left  to  the  right  and 
then  from  the  right  to  the  left  and  so  on.  The  data  acquisition  direction 
for  even-number  hybrid  scans  is  just  the  reverse  of  that  for  odd- 
number  hybrid  scans. 
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1st  hybrid  2nd  hybrid  3rd  hybrid 


RF  is  the  same  for  all  hybrid  scans. 

Figure  5.9.  The  pulse  sequence  foral28x  128  odd-number  hybrid  EPI.  The 
drawing  shows  the  sequence  for  the  three  hybrid  scans  used  in  the 
128  x  128  odd-number  hybrid  EPI.  The  time  delay  and  polarity  change 
of  the  X  gradient  pulses  are  shown. 
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Figure  5.11.  The  looping  technique  for  the  EPI  based  spin  tagging  sequence. 
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The  magnetization  transfer  pulse  is  a  lengthy  RF  pulse,  which  could  double  the  optr 
period  of  the  fast  SPGR  sequence  and  thus  destroy  the  time  constraint  of  the  spin 
tagging  sequence.  The  magnetization  transfer  RF  pulse  selectively  excites  the  spins  of 
the  macromolecules  such  as  protein  by  a  90°  flip  angle.  These  spins  have  a  Larmor 
frequency  offset  (either  lower  or  higher)  from  that  of  spins  in  water  due  to  the 
different  micro-environments  and  thus  can  be  selectively  excited.  Two  gradient 
spoiling  pulses  in  the  X  and  Y  directions  follow  this  RF  pulse  to  spoil  the  transverse 
magnetization  of  the  spins.  The  data  acquisition  RF  pulse  is  applied  immediately  after 
the  spoiling  pulses.  Only  the  spins  of  water  are  left  to  be  excited  by  this  RF  pulse 
because  longitudinal  magnetization  of  the  spins  from  macromolecules  has  not  yet 
recovered.  The  EPI  train  is  followed  as  in  a  regular  interleaved  EPI  sequence  for  data 
acquisition  (Figure  5.12).  The  amount  of  magnetization  transfer  measured  is  still  a 
relative  figure  rather  than  a  quantitative  one.  The  Larmor  frequency  of  the  spins  of  the 
macromolecules  can  be  randomly  distributed  over  a  large  range  around  the  Larmor 
frequency  of  the  spins  in  water.  The  ability  of  the  magnetization  transfer  RF  pulse  and 
its  associated  spoiling  gradient  pulses  are  also  limited  in  removing  the  effect  of  the 
magnetization  transfer  completely.  Due  to  hardware  problems,  this  pulse  sequence  has 
only  been  implemented  using  computer  simulation. 

This  chapter  has  provided  a  detailed  description  on  the  fast  SPGR  based  and 
the  EPI  based  arterial  spin  tagging  pulse  sequences  that  were  developed  and 
implemented  for  this  dissertation.  The  design  of  the  optimized  RF  inversion  and  data 
acquisition  pulses  has  also  been  discussed.  The  goal  is  to  use  these  pulse  sequences  to 
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Figure  5.12.  The  EPI  based  spin  tagging  sequence  with  MT  pulses  incorporated. 
The  odd-number  interleave  EPI  technique  discussed  is  implemented  in  data 
acquisition.  The  upper  figure  shows  the  timing  diagram  for  each  interleave  at 
acquisition.  A  180°  degree  RF  pulse  is  followed  by  N  repeated  EPI  data 
acquisitions  (as  shown  at  the  lower  figure). 
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capture  the  ideal  spin  recovery  after  the  application  of  the  spin  tagging  inversion  RF 
pulse  under  the  Selective  or  Non-selective  tagging  condition.  The  images  acquired  can 
be  used  to  extract  the  7Ys  under  the  Selective  and  Non-selective  conditions,  and  then 
the  perfusion  parameter / IX  can  be  calculated.  The  techniques  of  data  extraction  and 
processing  will  be  covered  in  the  next  chapter. 


Chapter  6 

Data  Processing  Techniques 


For  each  slice  location  to  be  analyzed  after  the  necessary  images  are  acquired 
using  the  arterial  spin  tagging  sequence,  the  pixel  values  across  all  dynamic  images 
are  used  to  calculate  the  T\  value  for  the  corresponding  pixel  location.  This  chapter 
will  focus  on  the  discussion  of  the  semi-log  linear  regression  method  used  to  extract 
the  T\  values  from  the  images.  Error-reduction  techniques  used  to  support  this  semi¬ 
log  linear  regression  method  will  also  be  discussed. 

6.1.  Semi-log  Linear  Regression  Method  for  T\  Calculation 

A  statistical  semi-log  linear  regression  technique  has  been  developed  based  on 
the  comparison  between  the  signal  from  undisturbed  spins  under  the  Regular 
condition,  and  the  signal  from  inverted  spins  under  the  Selective  or  Non-selective 

condition  (Figure  6.1).  The  b  value  (the  inverse  of  T\s  or  7i„)  is  equal  to  the  negative 

i 

slope  of  the  regression  line.  The  semi-log  linear  regression  equations,  which  will  be 
derived,  are  shown  below  based  on  two  arbitrary  data  points  (the  z'th  and y'th  data 
points)  of  the  T\  recovery  curve: 


b  =_L  =  hCflg  -  Ssei( 0)  -  HSng  -  SJM 
'  Tu  ’  TIbr(i)  -  TIbr(j) 


J_  =  ln(Sreg  -  Snon( Q)  -  In (Sng  -  Snon(j)) 
4  "  TIbr(i)  -  TIbrU) 


(6.1) 


(6.2) 
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The  steady-state  signal  intensity 
difference  in  log  scale  at  different  77 


* 


Figure  6.1.  An  example  of  semi-log  linear  regression  for  Tx  calculation.  This 
figure  shows  an  example  of  the  use  of  the  semi-log  linear  regression  method  in 
calculating  Tx  values  (TXs  at  Selective  condition,  and  r,„  at  Non-selective 
condition)  and  the  relative  perfusion  rate  / /A.  The  calculation  results  are  based  on 
the  number  of  points  left  after  the  application  of  the  minimum-mean-square-error 
technique.  TIbr  is  the  starting  time  of  image  acquisition  after  the  spins  are  inverted. 
Sreg  is  the  average  signal  intensity  at  one  pixel  from  the  last  three  images  at  the 
Regular  condition.  Sse,  is  the  signal  intensity  at  one  pixel  from  one  of  the  steady- 
state  images  at  the  Selective  condition,  S„on  is  the  signal  intensity  at  one  pixel 
from  one  of  the  steady-state  images  at  the  Non-selective  condition. 


where 


Sreg  =  the  average  signal  intensity  at  a  pixel  location  from  the  steady-state 
images  under  the  Regular  condition, 

TIbr{i)  or  TIbr(j)  =  the  starting  time  for  image  data  acquisition  after  spin 
inversion  at  the  rth  or y'th  timing  condition, 

SseAj)  or  Ssei(j)  =  the  signal  intensity  at  a  pixel  location  for  a  steady-state  image 
at  the  rth  or y'th  timing  condition  under  the  Selective  condition, 
Sn0„(i )  or  Sn0„(j ) =  the  signal  intensity  at  a  pixel  location  for  a  steady-state 
image  at  the  rth  or yth  timing  condition  under  the  Non- 
selective  condition. 

The  value  of f  IX  is  the  difference  between  the  inverse  of  the  T\  values  as  expressed  in 
Eq.  (4.7),  and  is  restated  here: 

/ 

J  =  bs  '  K  (4.7) 

i 

The  linear  regression  equations,  Eqs.  (6.1)  and  (6.2),  are  derived  below. 

In  this  derivation,  the  following  ideal  situation  is  assumed:  The  read-out  RF 
pulse  causes  such  a  small  flip  angle  that  its  continuous  application  for  all  the  phase¬ 
encoding  steps  (in  either  2D  fast  SPGR  or  EPI)  does  not  disturb  the  natural  recovery 
of  the  inverted  spins  under  the  tagging  conditions  and  the  equilibrium  states  of  the 
spins  under  the  non-tagging  condition. 
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Given  this  assumption,  under  a  spin  tagging  condition,  the  steady-state 
longitudinal  magnetization  at  the  end  of  each  TR  period  just  before  the  next  inversion 
pulse  can  be  derived  from  Eq.  (2.34)  with  flip  angle  a  =  180°.  The  result  is: 

=  K  (1  -  e-b  m )  -  e-‘' m  (6-3) 

where, 

=  the  steady-state  longitudinal  magnetization  at  the  end  of  each  TR 
period, 

Mq  =  the  longitudinal  magnetization  under  fully  relaxed  conditions. 


—  for  the  Selective  condition 
T 

6  =  j  j1 

—  for  the  Non  -  selective  condition 
TR  =  the  repetition  time. 

Thus, 


^ lOss  M0 


1-e 


-bTR 


l  +  e-*TR-  (6.4) 

The  longitudinal  magnetization  at  time  t  after  the  spin  inversion  at  steady-state  is 


M,(t)=M0(l-e-b')-Ml0sse-b' 

.  1  -  e'bTR 

=  M0(l-e-b')-M0— ^e- 

i  _  P~bTR 

=  M0(\-e~b'  -j—^e-b') 
2e~bl 

=  Mod  -T— 35F) 


or 
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(6.5) 


As  a  result,  the  signal  measured  at  echo  time  TE  is 


0p~bTI  JE 

S,  =  A/0(l-— ^-)sin(0)e  « 
l  +  e 


where, 


77  =  the  time  between  the  inversion  pulse  and  the  effective  center  of  k  space, 
TE  -  the  time  of  echo  at  each  read-out  period, 

6  =  the  read-out  RF  flip  angle. 

Under  the  spin  non-tagging  condition  (Regular  condition),  based  on  the  given 
assumption,  the  longitudinal  magnetization  Mre^t)  at  any  time  t  stays  at  equilibrium. 


Specifically, 


MnJt)=M0 


The  signal  measured  at  echo  time  TE  becomes 


Thus, 


Sng  =  M0  sin(0)e 


S,  =  2e~m 

C  1  t  -bTR 

Sreg  l  +  e 


Then, 


S,  =  2e~bn 
Sng  ~  \  +  e-b™ 
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s„g-s(_  2e-bn 

sng  i+e-ira 

Thus, 

S^-SXi)  _e~knw 

S^-S.U)  e~b71u) 

_  e-HJi(i)-nu)) 

where, 

77(i)  =  the  time  of  inversion  at  the  rth  77  cycle, 

SXO  =  the  signal  at  a  pixel  when  TI(i)  is  used, 

TI(j)  =  the  time  of  inversion  at  the  y'th  77  cycle, 

S/J)  =  the  signal  at  a  pixel  when  TI(j)  is  used. 

Then, 

In (Sng  -  S,  (/))  -  In (Sng  -  S,  (/))  =  -6(77(0  -  77(7)) 

Thus, 

^  ln(Sreg-^(0)~ln(^g-^(7)) 

r/(0-77(y) 


(6.10) 


(6.11) 


The  above  equation  shows  the  semi-log  linear  regression  relationship  between  (S^g  - 
St)  and  77.  It  can  be  applied  to  both  Selective  and  Non-selective  tagging  conditions. 
The  exact  locations  of  the  effective  77 s  are  uncertain,  but  they  should  contain  the 
same  amount  of  shift  from  their  apparent  positions  (the  times  from  the  spin  inversion 
to  the  starting  point  of  data  acquisition).  Because  of  the  subtraction  of  77(i)  and  TI(f), 
the  uncertainty  of  the  effective  TI  is  cancelled  out.  Therefore,  the  TI(i)  and  TI(j)  in  Eq. 
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(6.1 1)  can  be  replaced  with  ThAj)  and  TIbAj),  which  have  been  defined  in  Eqs.  (6.1) 
and  (6.2).  Then  Eq.  (6.1 1)  becomes 


ln(S„g  ~  S,  (0)  ~  bi(Sreg  -  S,  (j)) 
TIbr(i)-TIbrU ) 


(6.12) 


Eq.  (6.12)  can  be  expanded  to  a  higher  level  of  detail.  The  result  is  Eqs.  (6.1)  and  (6.2) 
that  have  been  stated  in  the  early  part  of  this  chapter,  and  can  be  restated  below: 


bs=  —  = 

Tu 


In -  Ssel(i ))  -  ln(Sreg  -  Ssel(j )) 


TIbr(i)  -  TIbr(j ) 


(6.1) 


b  = 


1 


ln (S^  -  Snon(i))  -  HSng  -  Snon(j)) 


”  T, 


TIbr(i)  -  TIbr(j) 


(6.2) 


Eq.  (6.1)  also  allows  for  the  elimination  of  the  effect  of  the  geometric 
mismatch  between  the  inversion  and  data  acquisition  slice  profiles  on  Tis  calculation. 

The  RF  excitation  pulses  used  for  the  spin  inversion  and  the  data  acquisition  are 

* 

different.  Under  the  Selective  condition,  geometric  mismatch  occurs  between  the 
inversion  and  data  acquisition  slice  profiles.  (Under  the  Non-selective  condition,  this 
is  not  a  problem  because  the  RF  inversion  pulse  inverts  spins  everywhere.)  Eq.  (6.1) 
can  be  used  to  address  this  problem  as  discussed  below: 

The  geometric  mismatch  between  the  inversion  and  read-out  profiles  can  lead 
to  signal  contamination  by  signal  from  non-inverted  spins.  Thus  the  apparent  signal 
from  the  Selective  condition  can  be  modeled  as 
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(6.13) 


=(1  -err)-S„,+err  -Sng 


where, 


S.,,  =  the  signal  of  the  same  pixel  actually  measured, 

lapp 

Ssei  =  the  signal  of  a  pixel  if  there  is  no  slice  profile  mismatch, 
err  =  the  fraction  of  error  due  to  slice  profile  mismatch. 


Then, 


Thus, 


-^*-  =  (1  -err)^-  +  err 

^ reg  ^reg 


And 


Ssel  _ 


-err 


sng  l~err 


I  ^sel  _  j  _  ^reg 


-err 


reg 


1  -err 


Sreg  ~  Ssel 


'reg 


1 

1-err - —  +  err 

_ Sreg 

1  -err 


1- 


Jse!„ 


reg 


1  -err 
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(6.14) 


j  x,<rp  '  ' 

Sng  ~Ssel(  0 _ Sftg 

sng-sse,(j)  ~  1  ^(y) 

_  Srrg-Ss'l^O) 

Srtg~Sselm,U) 


Comparing  with  Eq.  (6.10),  we  have 


(0  _  e-b,(Tl(i)-TIU)) 
S«g-Sseljj)- 


where, 


Thus, 


1  1  / 

Tu  Tx  +  X 


the  b  value  under  the  Selective  condition . 


b  ^ng-sselji))-HSng-sstljj)) 

s  77(0-770) 

HSKg-Ss^(i))-HSng-SselJj))  (6.15) 

'  TIbr(i)-TIbr(j) 

Eq.  (6.15)  corresponds  exactly  to  Eq.  (6.1).  In  other  words,  by  using  the  semi-log 
linear  regression  method,  the  problem  due  to  geometric  mismatch  between  the 
inversion  and  read-out  slice  profiles  is  eliminated.  For  simplicity  of  the  notations,  the 
subscript  "app"  has  not  been  used  in  Eqs.  (6.1)  and  (6.2),  but  the  Ssei  and  Snon  in  these 
two  equations  actually  represent  the  apparent  signals  at  a  pixel. 

6.2.  Error-Reduction  Techniques  for  T\  Calculation 

Error-reduction  techniques  that  are  based  on  three  different  approaches  have 
been  used  to  support  the  semi-log  linear  regression  method  for  T\  calculation.  These 


approaches  involve  the  elimination  of  the  uncertain  data  points,  pixel  averaging,  and 
minimum  mean-square-error  optimization. 

Elimination  of  uncertain  data  points  during  data  processing  is  considered  first. 
In  the  fast  SPGR  based  spin  tagging  pulse  sequence,  four  images  are  acquired  in  each 
of  the  seven  77  cycles.  The  last  three  images  have  achieved  steady  state  and  are 
extracted  for  the  linear  regression  calculation.  This  is  the  case  for  the  Selective  and 
Non-selective  conditions,  thus,  there  are  a  total  of  21  points  for  the  calculation  of  7js 
or  T\n.  Five  images  are  acquired  under  the  Regular  condition.  Similarly,  the  last  three 
images  have  achieved  steady  state  and  are  extracted  for  the  linear  regression 
calculation.  All  images  acquired  are  represented  as  magnitude  images,  but  the  actual 
signal  under  the  zero  line  of  the  T\  recovery  curve  is  negative  (Figure  4.3).  These 

J 

image  signals  must  be  negated  in  order  to  use  the  above  equations.  The  three  points 
(acquired  at  a  same  77)  that  are  located  the^closest  to  the  zero  value  are  usually 
eliminated  because  it  is  difficult  to  decipher  whether  these  points  are  above  or  below 
the  zero  line.  As  a  result,  a  total  of  1 8  points  contribute  to  the  calculation  of  7j. 

The  7j  and  fIX  calculations  based  on  this  pulse  sequence  are  very  sensitive  to 
noise,  especially  displacement  artifacts.  Thus  the  second  and  third  error-reduction 
techniques  are  the  noise-reduction  techniques  that  have  been  implemented  to  reduce 
this  problem.  For  the  second  technique,  the  resolution  of  each  image  is  reduced  from 
256  x  256  to  a  lower  resolution,  for  example,  64  x  64,  by  pixel  averaging.  This  helps 
to  reduce  white  noise  as  well  as  the  effect  of  possible  small  object  displacement 
among  the  61  images  being  analyzed. 
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For  the  third  error-reduction  technique,  an  approach  with  minimum-mean- 
square-error  optimization  is  used  for  the  calculation  of  T\  values.  A  maximum  of  five 
noisy  data  points  can  be  eliminated  to  maximize  T\  calculation  precision,  and  to  avoid 
inaccuracy  caused  by  over-elimination  of  data  points  (Figure  6.2).  The  algorithm  can 
be  outlined  as  following: 

1 .  Start  with  all  1 8  points. 

2.  Calculate  the  b  value  (1/  T\)  and  its  standard  deviation,  the  slope  and  the  Y 
axis-intercept  of  the  line. 

3.  If  the  number  of  points  used  for  calculation  is  equal  to  13,  or  the  standard 
deviation  of  the  b  value  is  greater  than  the  previous  one  calculated,  stop 
and  go  to  Step  7,  else  continue. 

4.  Calculate  the  distance  between  the  line  and  each  of  the  points. 

5.  Eliminate  the  point  with  the  maximum  distance. 

6.  Go  to  Step  2. 

7.  Output  the  b  value  and  its  standard  deviation. 

8.  End. 

The  algorithm  presented  above  is  used  to  ensure  that  a  data  point  containing 
excessive  noise  is  not  used  for  T\  calculation.  However,  there  is  a  limit  to  the  number 

of  points  that  can  be  eliminated,  since  over  elimination  can  lead  to  bias  calculation. 

* 

Here,  a  maximum  of  five  points  can  be  eliminated:  three  points  from  an  undesirable  77 
timing  condition  can  be  eliminated  and  two  noisy  points  from  other  77  timing 
conditions  can  also  be  eliminated. 

For  the  EPI  based  spin  tagging  pulse  sequence,  steady  state  is  not  a  concern  in 
the  T\  calculation  because  all  dynamic  data  is  acquired  after  a  single  RF  inversion 
pulse.  All  dynamic  pixel  signals  measured  after  the  inversion  pulse,  except  the  one  or 
two  points  closest  to  the  zero  line,  are  used  for  the  T\  calculation.  Therefore,  the  above 
approaches  can  also  be  used  but  in  this  case  will  include  more  data  points. 
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Figure  6 
find  the  < 


This  chapter  has  focused  on  the  discussion  of  the  semi-log  linear  regression 
method  used  to  extract  the  7j  values  from  the  images.  Error-reduction  approaches  (the 
elimination  of  the  uncertain  data  points,  pixel  averaging  and  minimum  mean-square- 
error)  have  been  used  to  support  this  semi-log  linear  regression  method.  The  semi-log 
linear  regression  method  has  been  evaluated  and  compared  to  the  traditional  zero¬ 
crossing  method,  which  determines  the  7j  value  based  on  the  zero-crossing  point  of 
the  T\  recovery  curve,  and  the  curve  fitting  technique  based  on  the  exponential  nature 
of  the  7j  recovery.  Neither  the  zero-crossing  method  nor  the  curve  fitting  technique 
has  been  able  to  provide  an  accurate  T\  calculation  due  the  uncertainty  of  the  effective 
77  values  and  the  RF  slice  profiles  geometric  mismatch  that  have  been  discussed 
earlier  in  this  chapter.  The  performance  of  the  semi-log  linear  regression  method, 
along  with  the  three  error-reduction  approaches,  is  clearly  superior  to  the  zero¬ 
crossing  and  curve  fitting  techniques.  Thus  the  semi-log  linear  regression  method  and 

the  three  error-reduction  approaches.are  applied  to  all  the  data  used  in  this  work. 

* 

In  the  following  chapter,  an  image  processing  technique  that  was  used  to 
further  improve  the  7j  calculation  accuracy  through  object  displacement  artifact 
reduction  will  be  discussed.  Additionally,  image  visualization  and  analysis  techniques 
that  allow  for  the  application  of  the  calculation  results,  which  obtained  from  the  data 
processing  techniques  presented  in  this  chapter,  to  breast  cancer  detection  will  be 
presented. 
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Chapter  7 


Image  Processing  Techniques 

Motion  artifacts  are  always  problematic  for  dynamic  images  collected  over  a 
period  of  time.  How  should  they  be  corrected?  After  the  T\s,  7jn,  and  /// 1  values  have 
been  calculated  for  the  slices  being  investigated,  how  can  they  be  visualized  in  a  clear 
manner?  How  can  these  values  be  used  to  aid  clinicians  in  making  better  decisions  in 
identifying  cancerous  regions?  In  this  chapter,  techniques  developed  to  solve  these 
problems  are  discussed. 

7.1.  Binary  Segmentation  Technique  for  Displacement  Artifact  Correction 
An  important  image  processing  technique  that  has  been  developed  in  this 

1 

project  reduces  the  displacement  artifacts  that  can  cause  inaccuracies  in  7j 
calculations.  Both  breast  objects  (referring  to  the  two  objects  seen  on  each  individual 
image)  are  expected  to  be  in  the  same  relative  position  for  all  images  acquired  at  one 
slice  location.  Unfortunately,  the  relative  position  of  the  breast  objects  might  change 
during  the  three-minute  scan  in  the  fast  SPGR  based  spin  tagging  pulse  sequence. 
Researchers  would  like  to  know  the  amount  of  displacement  and  correct  it  when 
necessary.  A  typical  breast  perfusion  study  generates  61  images  at  each  slice  location, 
and  10  locations  are  analyzed,  giving  a  total  of  610  images.  Due  to  the  large  number 
of  images,  an  algorithm  requiring  intensive  computation  to  estimate  and  correct  object 
displacement  is  not  appropriate.  A  fast  and  effective  technique  is  needed.  This  goal 
has  been  achieved  by  estimating  dislocation  based  on  binary  segmentation  of  the  61 
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images  acquired  from  a  typical  breast  study  using  the  fast  SPGR  based  spin  tagging 
sequence.  When  the  EPI  based  spin  tagging  sequence  is  used,  the  images  are  collected 
in  much  less  time,  and  motion  artifacts  become  less  problematic.  The  displacement 
artifact  correction  technique  presented  here  can  still  be  applied  to  these  images  as  well 
as  be  generalized  to  other  types  of  dynamic  studies. 

Knowledge  about  a  dynamic  study  usually  can  help  to  improve  the  accuracy 
and  efficiency  of  image  processing,  and  certainly  reduces  the  complexity  of  the 
algorithm  design.  During  scanning,  the  subject  is  restricted  on  the  scanner  table,  and 
thus  rotation  of  the  breasts  is  virtually  impossible.  Only  lateral  displacement  in  the  X 
and  Y  directions  needs  to  be  considered.  Fortunately  the  displacement  is  not  large, 
most  likely  within  a  range  of  +/-1  cm.  The  approximate  shapes,  sizes  and  positions  of 
the  objects  are  known.  Although  signal  intensity  varies  greatly  from  one  image  to 
another,  with  some  images  having  poor  object-background  separation,  the  breast 

objects  are  still  separable  from  the  background.  This  knowledge  leads  to  the 

* 

implementation  of  a  binary  segmentation  technique  discussed  below. 

The  61st  image  is  obtained  at  the  steady  state  of  the  Regular  condition,  and  is 
proton-density  weighted.  It  has  a  good  object-background  separation  and  thus  binary 
segmentation  is  clean,  so  it  is  used  as  the  reference.  A  threshold  is  easily  chosen  for 
the  reference  image  to  separate  the  objects  from  their  background.  Based  on  this 
threshold,  almost  all  the  pixels  on  the  breast  objects  are  labeled  as  ones  and  almost  all 
the  pixels  on  the  background  are  labeled  as  zeros.  In  this  way,  a  good  reference 
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segmentation  image  is  created.  For  all  other  60  images,  once  a  threshold  is  determined 
as  discussed  in  the  next  paragraph,  binary  segmentation  images  are  created  similarly. 

Since  the  reference  image  has  a  very  clean  object-background  separation,  the 
total  number  of  pixels  in  the  two  breast  objects  (population  of  object  pixels)  is 
approximately  equal  to  the  mass  of  the  binary  segmentation  image  (or  the  total  number 
of  ones).  When  one  of  the  other  60  images  is  analyzed,  a  histogram  of  pixel  values  is 
created  first.  The  pixel  value  that  separates  the  populations  of  object  and  background 
pixels  is  the  threshold  used  for  binary  segmentation  [51]. 

Cross  correlation  is  done  between  the  segmentation  image  to  be  realigned  and 
the  reference  segmentation  image.  Maximum  cross  correlation  occurs  when  the  two 
images  overlap.  Since  the  range  of  breast  object  displacement  is  short,  only  a  local 

i 

cross  correlation  is  needed  (Figure  7.1).  This  greatly  reduces  the  computation  time. 
This  also  takes  advantage  of  the  fact  that  the  pixel  values  at  the  outer  parts  of  the 
objects  can  be  separated  better  from -their  backgrounds  than  those  at  the  inner  parts  of 
the  objects. 

To  further  speed  up  the  computation,  two  iterations  are  used  to  compute  the 
cross  correlation: 

Range  of  computation:  -10  to  10  pixels  in  both  X  and  Y  directions, 

Iteration  #1:  -8,  -4,  0, 4,  8, 

Iteration  #  2:  -2,  -1,0, 1,2. 

The  displacement  estimation  and  correction  algorithm  was  implemented  using 
the  C  language  and  was  evaluated  on  an  SGI  02  (195  MHz,  MIPS,  R10000) 
workstation.  For  the  analysis  of  the  61  images,  the  run  time  was  2.2  minutes. 
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Reference  image 

Image  to  be  realigned 
Region  for  local  cross  correlation 


Figure  7.1.  Object  displacement  estimation  based  on  local  cross 
correlation.  Local  cross  correlation  is  performed  between  the  binary  image 
for  realignment  against  the  reference  binary  image.  The  maximum 
correlation  occurs  when  two  objects  overlap. 


Figure  7.2  shows  the  displacement  correction  applied  to  an  artificially  introduced 
displacement.  The  average  mis-calculation  of  object  displacement  was  0.39  pixel.  This 
result  was  found  to  be  superior  to  the  edge-detection  technique  and  other  binary 
segmentation  techniques  that  had  been  attempted  by  the  author.  Since  the  threshold  is 
automatically  determined,  this  technique  has  high  portability. 

7.2  Image  Visualization  and  Segmentation 

After  the  necessary  displacement  correction,  the  values  of  7js,  7j„  and / // 1  at 
different  pixel  locations  on  each  slice  can  be  calculated  based  on  the  semi-log  linear 
regression  technique.  Images  can  be  formed  based  on  these  values.  They  are  called 
feature  images  because  they  represent  the  characteristics  of  breast  tissue.  The  most 
useful  feature  images  for  tissue  characterization  are  the  7j  image  based  on  the  Tin 
values  in  the  units  of  sec,  the  perfusion  image  based  on  the  / IX  values  in  the  units  of 
sec'1  and  the  error  perfusion  image  based  oh  the  standard  deviations  of  the / IX  values 
also  in  the  units  of  sec'1.  These  feature  images  contain  pixels  with  a  wide  range  of 
values.  But  the  values  of  7j,/  IX  and  the  standard  deviation  of f  IX  of  breast  tissues 
should  be  within  reasonable  ranges.  Pixel  values  below  or  above  certain  values  are  not 
useful  and  can  be  misleading.  To  display  these  feature  images,  an  appropriate  positive 
shift,  scaling  and  window  size  should  be  applied  (Figure  7.3)  [52].  Linear  scaling  is 
applied  to  preserve  the  characteristics  of  the  data.  The  details  of  these  manipulations 
are  shown  in  Table  7.1  below: 
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(e) 


Figure  7.2.  Displacement  restoration  based  on  binary  segmentation,  (a)  Reference 
(Image  61).  (b)  The  segmented  binary  image  of  Figure  (a)  with  threshold  100.  (c) 
Displacement  was  artificially  introduced  to  Image  8.  The  left  object  was  shifted  to 
the  right  by  four  pixels  and  shifted  up  by  five  pixels.  The  right  object  was  shifted  to 
the  left  by  six  pixels  and  shifted  down  by  six  pixels,  (d)  The  segmented  binary 
image  of  Figure  (c)  with  threshold  70.  (e)  The  restored  image  of  Figure  (c). 
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Figure  7.3.  Feature  image  display  conversion.  Before  a  feature  image  is 
displayed,  it  needs  to  be  processed  to  provide  a  meaningful  display  range 
and  a  reasonable  contrast  level  by  contrast  stretching  with  appropriate 
cutoff,  offset  shift  and  scaling. 
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Table  7.1.  The  Contrast  Stretching  Parameters  Used  for  the 
Display  of  the  Three  Feature  Images 


Feature  images 

T\  image  (sec) 

/ IX  image  (sec'1) 

STD  of  fIX  image  (sec'1) 

Possible  range  of 
values 

from  0  to  4 

from  -4  to  4 

from  -5  to  5 

Range  of  interest 

from  0  to  1.5 

from  -0.5  to  0.8 

from  -0.5  to  0.5 

Positive  shift 

0 

0.5 

6.5 

Upper  cutoff 

1.5 

0.8 

0.5 

Lower  cutoff 

0 

-0.5 

-0.5 

Scaling  factor 

1000 

1000 

1000 

A  visualization  computer  program  called  BreastView  has  been  developed  by 
the  author  on  an  SGI  workstation  based  on  C  language,  X  Window  System  and  Motif 
1.2.  This  program  will  be  used  to  interpret  these  three  feature  images  for  the 
recognition  of  suspicious  cancer  regions.  Since  a  sizable  tumor  extends  to  several 
pixels,  these  pixels  should  have  similar  tissue  characteristics.  Specifically,  they  should 
have  high  T\  values  that  are  approximately  equal  to  each  other,  and  high  perfusion 
rates.  When  a  pixel  is  found  to  have  a  high  T\  and  a  high / IX,  and  the  standard 
deviation  of  the  f/X  is  low,  its  surrounding  pixels  within  the  same  slice  and  across 
slices  would  be  checked  to  see  if  they  have  similar  characteristics.  The  level  of 
suspicion  for  this  pixel  is  calculated  based  on  the  percentage  of  its  surrounding  pixels 
(including  this  central  pixel)  containing  the  similar  suspicious  characteristics.  For  a 
pixel  on  an  edge  slice,  there  is  a  total  of  1 8  pixels  to  check;  otherwise,  there  is  a  total 
of  27  pixels  to  check.  By  applying  this  method  to  all  pixels  in  all  slices,  maps  of 
suspicion  are  created.  If  the  pixel  has  a  total  level  of  suspicion  exceeding  a  cutoff 
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value,  such  as  20%,  this  pixel  would  be  identified  as  a  cancer  pixel.  The  goal  is  to 

develop  an  analytical  tool  for  breast  cancer  detection.  The  exact  cutoff  will  depend  on 

the  judgement  of  researchers  and  clinicians.  The  suspicion  level  cutoff,  as  well  as  the 

cutoffs  for  f/X,  T\  and  the  standard  deviation  of f  IX  should  be  obtained  based  on  a 

large  set  of  studies.  Based  on  the  small  set  of  studies  that  have  been  done,  the 

approximate  cutoffs  that  can  be  used  to  identify  cancer  pixels  are  listed: 

7j:  >  0.5  sec 

f/X :  >  0.1  sec'1 

STD  of//A:  <  0.1  sec'1. 

Suspicion  level  threshold:  >  20  %. 

Researchers  will  be  able  to  use  the  BreastView  program  to  manipulate  the 
above  cutoffs  and  to  observe  the  change  of  the  suspicion  maps.  This  visualization 
program  can  simultaneously  display  the  feature  images  and  clinical  images. 

Parameters  T\,f IX,  and  the  standard  deviation  of f  IX  at  any  pixel  location  of  the  breast 
objects  can  be  found  by  a  simple  click  of  the  mouse.  If  a  dynamic  contrast 
enhancement  study  fias  been  done,  a  click  of  the  mouse  can  also  activate  the  program 
to  plot  the  dynamic  change  of  signal  intensity  at  this  location.  The  suspicious  cancer 
pixels  can  also  be  mapped  to  high-resolution  traditional  clinical  images  (256  x  256), 
which  are  always  acquired  for  comparison  during  each  subject  study.  This  mapping 
provides  a  radiologist  with  the  traditional  view  of  the  breast  images  with  the 
suspicious  cancer  regions  identified.  BreastView  is  designed  to  allow  easy  pixel  cross¬ 
registration  and  comparison  between  results  from  different  MR  imaging  techniques. 
This  would  provide  rich  information  for  clinicians  in  making  more  accurate  decisions. 
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Appendix  B  provides  more  detail  information  on  the  installation  and  operation  of  this 
program. 

Images  from  a  patient  with  a  suspicious  mass  found  by  palpation  in  her  right 
breast  were  used  to  verily  the  image  processing  methods  presented.  Figure  7.4  shows 
the  three  feature  images  created  along  with  an  image  acquired  under  the  Regular 
condition  of  the  breast  sequence  for  comparison.  All  three  feature  images  have  the 
resolution  of  64  x  64.  The  T\  image  has  a  good  SNR.  It  clearly  defines  the  distribution 
of  T\  values  on  the  slice.  The  brighter  areas  at  the  inner  regions  of  the  breasts  indicate 
higher  T\  values  and  thus  higher  water  content.  This  perfusion  image  appears  to  have  a 
good  SNR  though  sometimes  the  perfusion  image  is  quite  noisy.  This  happens  because 
the  intensity  of  the  white  noise  is  within  the  range  of  the / IX  values.  A  pixel  would  not 
stand  out  unless  it  had  a  very  high  f  IX  value.  For  both  the  T\  and  the  perfusion  images, 
some  pixels  at  the  edge  of  the  breasts  have  a  high  intensity,  but  they  are  just  edge 
artifacts  due  to  high  calculation  error  at  the  edges.  Since  pixel  values  at  the  edges  are 
sensitive  to  breast  object  displacement  as  well  as  ringing  artifacts,  high  calculation 
error  can  happen  at  these  pixels.  The  error  levels  of  the / IX  calculation  can  be  found 
on  the  error  perfusion  image.  Pixels  with  low  T\  values  generally  contain  more  error  in 
the  calculation  of  perfusion  because  the  equations  (Eqs.  (6.1)  and  (6.2))  used  for  the 
calculation  of  perfusion  are  more  sensitive  to  the  fluctuation  of  T\  at  low  T\  values. 
Fortunately,  this  does  not  cause  much  trouble  because  cancer  tissue  should  have 
relatively  high  T\  values  (>  500  ms). 
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Figure  7.4.  Three  feature  images  derived  for  one  slice  location  in  one 
subject  study,  (a)  The  last  image  under  the  Regular  condition  is 
displayed  for  comparison.  The  three  feature  images  are  (b)  the  7j  image, 
(c)  the  perfusion  image  based  on / IX,  and  (d)  the  error  perfusion  image 
based  on  the  standard  deviation  of f  IX. 
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Figure  7.5  shows  the  mapping  of  suspicious  cancer  pixels  onto  high-resolution 
traditional  clinical  images  according  to  the  methods  discussed.  The  suspicious  regions 
identified  by  the  pure  white  pixels  on  the  right  breasts  in  Slices  6,  7  are  considered  one 
connected  region.  These  pure  white  pixels  are  actually  red  in  the  visualization  program 
BreastView.  This  suspicious  region  is  at,  or  right  next  to,  the  benign  mass  that  appears 
bright  in  both  clinical  images.  The  malignancy  of  this  suspicious  region  was  later 
confirmed  by  needle  biopsy.  In  general,  for  positive  identification  of  cancer,  these 
suspicious  regions  identified  in  the  breast  sequence  should  be  compared  with 
traditional  clinical  MR  images,  as  well  as  X-ray  mammography  and  ultrasound 
images.  It  is  often  unclear  whether  abnormal  masses  found  by  palpation,  traditional 
MRI  and  X-ray  mammography  are  malignant  or  not.  Suspicious  cancer  regions 

i 

identified  by  this  method  provide  important  information  in  confirming  the  malignancy 
of  the  abnormal  masses. 

7.3.  Bayesian  Technique  and  ROC 

* 

The  suspicious  cancer  pixel  identification  based  on  thresholds  provides  the 
basic  approach  for  cancer  recognition  in  this  work.  The  number  of  cases  studied  here 
is  limited.  However,  if  a  large  pool  of  case  studies  is  available,  the  Bayesian  decision 
technique  [53]  in  cancer  identification  can  be  applied.  The  Bayesian  technique  to  be 
discussed  can  be  extended  to  the  ROC  (receiver  operating  characteristic)  curve  that  is 
commonly  used  for  the  evaluation  of  a  diagnostic  algorithm  in  radiology. 
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Proton  Density  Weighted,  Slice  6 


T2  Weighted,  Slice  6 


Proton  Density  Weighted,  Slice  8  T2  Weighted,  Slice  8 


Figure  7.5.  Suspicious  cancer  region  color  mapping.  The  identified 
suspicious  cancer  pixels  are  mapped  onto  the  high-resolution  traditional 
clinical  MR  images  according  to  their  corresponding  slice  locations.  These 
images  are  from  three  consecutive  3mm  slices  with  inter-slice  distance  of 
10  mm.  The  pure  white  mapping  pixels  are  actually  red  in  the  visualization 
program  BreastView. 


This  Bayesian  technique  can  be  outlined  as  follows: 

Let  vector  space  x  =  (T\,f  IX,  STD  of f  IX)  for  a  ROI  (region  of  interest)  in  a 
breast.  The  risks  associated  with  a  given  vector  x  are 

Risk(Label  Benign  |  x)  =  r(Label  Benign  |  Malignant)  P(Malignant  |  x) 

+  r(Label  Benign  j  Benign)  P(Benign  |  x)  (7.1) 

Risk(Label  Malignant  |  x)  =  r(Label  Malignant  |  Malignant)  P(Malignant  |  x) 

+  r(Label  Malignant  |  Benign)  P(Benign  |  x)  (7.2) 

where, 

Risk(Label  Benign  |  x)  =  the  risk  value  to  label  a  ROI  "Benign"  for  a  given  x, 
Risk(Label  Malignant  |  x)  =  the  risk  value  to  label  a  ROI  "Malignant"  for  a 

given  x, 

r(Label  Benign  |  Malignant)  =  the  risk  factor  to  label  a  ROI  "Benign"  but 

actually  be  malignant, 

r(Label  Benign  |  Benign)  =  die  risk  factor  to  label  a  ROI  "Benign"  and  actually 

be  benign, 

r(Label  Malignant  |  Malignant)  =  the  risk  factor  to  label  a  ROI  "Malignant" 

and  actually  be  malignant, 

r(Label  Malignant  |  Benign)  =  the  risk  factor  to  label  a  ROI  "Malignant"  but 

actually  be  benign, 

P(Malignant  |  x)  =  the  probability  of  being  malignant  for  a  given  x, 

P(Benign  |  x)  =  the  probability  of  being  benign  for  a  given  x. 

Generally  there  is  not  any  risk  for  correct  labeling,  and  thus 
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r(Label  Benign  |  Benign)  =  0 


and 

r(Label  Malignant  |  Malignant)  =  0 
Then  Eq.  (7.1)  and  (7.2)  become 

Risk(Label  Benign  |  x)  =  r(Label  Benign  |  Malignant)  P(Malignant  |  x)  (7.3) 

Risk(Label  Malignant  |  x)  =  r(Label  Malignant  |  Benign)  P(Benign  |  x)  (7.4) 


Since, 


P(Malignant|x) 


p(x|  Malignant)P(Malignant) 
P(x) 


(7.5) 


and 


p^ni^^^PWBemg^Ben^ 

p(x) 


(7.6) 


where, 

* 

p(x  |  Malignant)  =  the  probability  density  of  vector  x  given  that  the  ROI  is 
malignant, 

P(Malignant)  =  the  probability  of  being  malignant  for  a  ROI 

p(x  |  Benign)  =  the  probability  density  of  vector  x  given  that  the  ROI  is  benign, 

P(Benign)  =  the  probability  of  being  benign  for  a  ROI, 

p(x)  =  the  probability  density  of  vector  x. 
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Risk  Ratio 


Risk(Label  Benign  |  x) 

Risk(Label  Malignant  |x) 
r(Label  Benign  |  Malignant)  P(Malignant  |  x) 
r(Label  Malignant  |  Benign)  P (Benign  |  x) 
r(Label  Benign  |  Malignant)  p(x  |  Malignant)  P(Malignant) 
r(Label  Malignant  j  Benign)  p(x  |  Benign)  P(Benign) 


(7.7) 


The  ratio  ^en*gn  I  Malignant)  .g  cjj0sen  ^g  clinician  or  researcher. 

r(Label  Malignant  |  Benign) 

The  ratio  1  ^nant-  is  composed  from  a  training  data  set.  The  ratio 


p(x  |  Benign) 


P(Malignant) 


can  be  obtained  from  literature  or  estimated  from  the  training  data  set; 


P(Benign) 

but  the  source  of  this  ratio  obtained  must  correspond  with  the  source  of  subjects  used 


in  the  breast  study.  Estimation  of  the  two  ratios,  ^-X-  ^  --gnan^  and  } , 

p(x  |  Benign)  P(Benign) 

requires  a  large  training  data  set  to  be  statistically  significant. 

If  the  Risk  Ratio  is  greater  than  1,  the  ROI  would  be  labeled  "Malignant". 

Otherwise,  the  ROI  would  be  labeled  "Benign".  Or 


Risk  Ratio 

ROI  Label 

>  1 

Malignant 

<  1 

Benign 

The  values  of  risk  factors  are  determined  by  the  clinician  or  the  researcher 
based  on  his/her  experience.  For  example,  if  the  clinician  thinks  that  it  is  too  risky  to 
label  a  malignant  mass  to  be  "Benign",  he/she  would  give  a  high  value  for  r(Label 
Benign  |  Malignant),  and  a  relative  low  value  for  r(Label  Malignant  |  Benign). 


L 
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Choosing  the  different  values  of  the  risk  factor  ratio  1  Malignant) 

r(Label  Malignant  j  Benign) 

equivalent  to  choosing  different  values  of  the  decision  threshold.  The  pairs  of  false 
positive  and  true  positive  fractions  resulting  from  the  different  values  of  the  risk  factor 
ratio  would  allow  for  the  creation  of  an  ROC  curve,  which  can  be  used  to  evaluate  the 
performance  of  the  technique  presented  [1]. 

This  chapter  has  covered  the  binary  segmentation  technique  used  for 
displacement  artifact  correction,  the  methods  of  image  visualization  and  segmentation 
used  for  cancer  identification,  the  Bayesian  decision  technique  used  for  cancer 
recognition,  and  the  ROC  curve  used  for  statistical  analysis.  The  binary  segmentation 
technique  has  been  shown  to  be  effective  in  estimating  and  correcting  displacement 
artifacts.  The  availability  of  the  user-friendly  BreastView  program  that  was  written  to 
implement  the  image  visualization  and  segmentation  methods  presented  in  this  chapter 
to  clinicians  and  researchers  should  allow  for  more  confidence  in  the  identification  of 
cancerous  breast  tissues.  The  Bayesian  technique  can  be  used  for  the  identification  of 
cancerous  breast  tissue  if  a  large  pool  of  case  studies  is  available. 
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Chapter  8 

Results  and  Discussion 


The  performance  of  the  methods  based  on  the  arterial  spin  tagging  sequences 
to  measure  T\  and  perfusion  has  been  verified  with  a  flow  phantom,  meat  phantoms, 
normal  breasts,  and  breasts  with  abnormal  masses.  Most  of  these  studies  were 
performed  on  a  Signa  Advantage  1 .5T  GE  MR  system  and  the  rest  were  performed  on 
a  Signa  Horizon  LX  1.5T  GE  MR  system  (GE  Medical  Systems  Inc.,  Milwaukee).  A 
dual  phased  array  breast  imaging  RF  receiving  coil  (Medical  Advance,  Inc., 
Milwaukee)  was  used  for  the  breast  studies.  All  spin  tagging  studies  were  done  in 
coronal  view  with  the  subject  lying  down  in  a  prone  position.  The  coronal  view  is 

i 

parallel  to  the  body  and  is  expected  to  be  perpendicular  to  the  direction  of  maximum 
net  blood  flow  in  the  breasts. 

8.1.  Flow  Phantom  Study 

Perfusion  studies  using  a  controllable  flow  phantom  were  the  first  set  of  and 
the  most  important  verification  studies  done.  A  flow  phantom  was  assembled  from  a 
kidney  dialysis  tube.  This  dialysis  tube  is  composed  of  artificial  capillaries  that  allow 
the  exchange  of  molecules  during  dialysis  for  a  kidney  patient.  The  kidney  dialysis 
tube  is  a  very  good  tissue  model  for  perfusion  studies.  It  has  a  length  of  about  25  cm, 
an  internal  diameter  of  about  5  cm,  about  1000  capillaries  and  an  approximate  A,  of 
0.15  (the  fraction  of  the  "blood"  volume  relative  to  the  total  volume  including  the 
"blood"  and  the  "interstitial  fluid").  This  "blood"  volume  is  assumed  to  be  the  total 
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volume  flowing  out  of  the  dialysis  tube  after  it  is  filled  up.  It  was  measured  to  be  20 
ml.  The  "interstitial  fluid"  is  assumed  to  be  the  total  volume  that  can  fill  up  the  space 
inside  the  tube  but  outside  of  the  capillaries.  It  was  measured  to  be  1 13  ml.  For  the 
convenience  of  physical  measurement,  the  perfusion  of  the  flow  phantom  was  only 
verified  in  one  direction,  instead  of  multiple  directions  that  occur  in  actual  tissue.  The 
flow  is  fed  through  long  rubber  tubes,  and  is  driven  by  an  external  electrical 
mechanical  pump  outside  the  magnet  room  to  avoid  hazards  associated  with  the  strong 
magnetic  field  (Figure  8.1).  The  fluid  used  in  the  phantom  was  purified  drinking  water 
doped  with  manganese  chloride  to  have  a  T\  value  of  approximately  895  ms,  which  is 
within  the  range  of  T\  value  of  body  fluid.  The  flow  rate  was  controlled  by  adjusting 
the  mechanical  pump  and  was  physically  measured  using  a  stop  watch  and  a  graduated 
cylinder. 

The  fast  SPGR  based  pulse  sequence  developed  has  been  optimized  for 

perfusion  measurements  by  compromising  between  scan  time  and  the  precision  of f  IP 

\ 

value.  All  phantom  studies  showed  that  there  was  a  linear  relationship  between f IP, 
and  the  actual  flow  applied  within  the  range  of  tissue  perfusion  (Figure  8.2).  In  these 
studies,  the  7jn  values  during  non-selective  inversion  at  four  different  flow  rates  were 
about  the  same  as  expected:  1.232  +/-  0.034  sec.  To  improve  measurement  sensitivity, 
the  "interstitial  fluid"  was  not  filled.  So  the  actual  value  of  the  P  in  this  phantom 
should  be  within  the  range  of  [0.15,  1].  It  can  be  estimated  based  on  the  physically 
measured  flow  rate  and  the / IP  value  measured  by  the  spin  tagging  technique.  For 
example,  at  a  physical  average  input  flow  rate  of  0.0222  ml/sec/cm2  through  the  cross 
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25  cm 


Figure  8.1.  The  setup  of  the  kidney  dialysis  phantom.  The  figure  shows 
the  mechanical  setup  of  the  kidney  dialysis  phantom  (1000  capillary 
tubes)  and  accessories.  The  flow  rate  can  be  controlled  through  the  pump. 
The  flow  rate  can  be  measured  through  the  fluid  receiving  end  of  the  tube 
after  it  is  pulled  out  from  the  fluid  container. 
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Dialysis  phantom  study  at  four  flow  rates 
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Figure  8.2.  Results  from  the  kidney  dialysis  phantom  perfusion  study. 
The  arterial  spin-tagging  sequence  was  applied  at  four  different  low 
flow  rates  with  0,  0.0222,  0.0640  and  0.0787  ml/sec/cm2.  The 
corresponding/ IX  values  calculated  were  0.0251,  0.138,  0.331  and 
0.409  sec'1. 


section  of  the  phantom,  the  /  IX  value  was  measured  to  be  0. 1 38  sec'1  through  the  spin 
tagging  technique.  Since  a  slice  thickness  of  0.3  cm  was  used  in  the  scan,  the  physical 
input  flow  rate  of  0.0222  ml/sec/cm2  is  equivalent  to  a  perfusion  value / of  0.0740 
sec'1,  which  is  obtained  by  the  simple  relationship  below: 

j.  flow  rate  (ml  /  sec  /  cm2) 
slice  thickness  (cm) 

Then  the  X  value  is  calculated  to  be  0.536  based  on  these / and / IX  values.  It  is  within 
the  range  of  [0.15, 1]  as  expected. 

8.2.  Meat/Fat  Phantom  Study 

The  meat/fat  phantom  study  is  used  to  verify  the  applicability  of  the  sequence 
to  tissue  similar  to  breast.  The  two  meat/fat  phantoms  are  made  of  pork  meat  enclosed 
with  both  pork  fat  and  chicken  fat  in  two  plastic  cups.  Each  plastic  cup  has  a  height  of 
6  cm,  a  top  diameter  of  9  cm  and  a  bottom  diameter  of  7  cm.  Since  the  tissue  in  this 
phantom  is  not  alive,,  it  should  have  nearly  zero  perfusion.  The  perfusion  measured 
should  just  be  caused  by  the  random  motion  of  water  molecules  within  the  tissue. 
Since  the  water  molecules  are  constrained  by  the  meat  and  fat  tissues,  their  motion 
should  be  very  limited.  Four  slice  locations  are  analyzed.  Thirty-eight  regions  of 
interest  (ROIs)  at  the  meat  and  21  ROIs  at  the  fat  from  these  four  slice  locations  are 
chosen  for  analysis  (Figure  8.3).  The  size  of  each  ROI  was  35  mm2.  Figure  8.4 
graphically  shows  the  distribution  of  7jn  and  / IX  in  meat  and  fat.  The  distribution  of 
Tin  and / IX  is  summarized  in  Table  8.1  below: 
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(a) 


(c) 


Figure  8.3.  The  meat/fat  phantom  images  obtained  with  the  fast  SPGR  based 
spin  tagging  sequence.  The  above  figures  show  images  obtained  at  (a) 
Selective,  (b)  Non-Selective  and  (c)  Regular  conditions  at  one  slice  location. 
All  images  are  displayed  using  the  window  level  of  Min/Max:  0/516.  The 
center  darker  area  in  Images  (a)  and  (b)  or  the  equivalent  brighter  area  in 
Image  (c)  is  meat  and  the  rest  of  the  area  in  the  two  objects  is  fat. 
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Table  8.1.  Meat-Fat  Characteristic  Comparison  in 
Meat/Fat  Phantom 


Mean  of 

Tin  (sec) 

STD  of 

Tin  (sec) 

Mean  of 
f/X  (sec1) 

STD  of 
/  IX  (sec*1) 

Fat 

0.363 

0.061 

-0.010 

0.116 

Meat 

0.598 

0.043 

0.022 

0.086 

The  Tin  of  meat  is  higher  than  that  of  fat.  The / IX  value  of  fat  is  closer  to  zero  than 
that  of  meat  because  the  water  content  in  fat  is  lower  than  that  in  meat.  The  higher 
standard  deviation  in  the  calculation  of  7j„  and / IX  is  due  to  the  higher  b  values  (the 
inverse  of  T\  values). 

8.3.  Human  Subject  Studies 

Eighteen  subjects  participated  in  the  breast  studies.  One  subject  also 
participated  in  the  first-pass  contrast  enhanced  dynamic  comparison  study.  All 
volunteers  signed  the  consent  form  of  research  agreement  that  had  been  approved  by 

V 

the  human  subjects  research  committee  at  University  of  California,  Davis.  At  each 
study,  some  simple  but  effective  procedures  were  used  to  reduce  motion  artifacts:  (1) 
the  volunteer  would  be  strapped  on  the  scan  table  in  a  prone  position  (Figure  8.5);  (2) 
the  volunteer  was  instructed  not  to  move  and  to  breathe  smoothly  during  scanning. 
Besides  the  examination  using  the  arterial  spin  tagging  pulse  sequence,  other 
examinations  using  traditional  clinical  techniques  were  also  done.  The  following 
procedure  was  followed  for  each  study: 
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Figure  8.5.  Subject  preparation  for  breast  imaging. 


1 .  2D  spin  echo  scan  at  the  axial  plane:  this  T\  weighted  multi-slice  scan  acts 
as  a  localizer  to  provide  the  overview  of  the  whole  breasts  and  neighboring 
tissues  and  helps  to  detect  the  existence  of  any  abnormal  mass. 

2.  2D  fast  spin  echo  scan  at  the  coronal  plane:  this  proton  density-weighted 
and  72-weighted  scan  provides  the  traditional  clinical  breast  images.  The 
resulting  images  are  used  to  compare  with  the  images  obtained  from  the 
arterial  spin  tagging  sequence. 

3.  Arterial  spin  tagging  pulse  sequence  at  the  coronal  plane:  about  10  slices 
centering  around  the  abnormal  masses  that  have  been  identified  by  the 
patient  and/or  the  patient's  physician  or  have  been  identified  in  the  two 
previous  clinical  scans.  For  normal  subjects,  about  10  slices  locations  are 
chosen  at  the  central  regions  of  the  breasts  with  inter-slice  distance  of  10 
mm.  The  10  slices  should  be  at  the  same  locations  as  the  ones  obtained  in 
the  2D  fast  spin  echo  scan,  and  should  have  the  same  field  of  view.  This 
would  allow  for  direct  comparison  between  clinical  techniques  and  the 
arterial  spin  tagging  technique. 

If  the  first-pass  contrast  enhanced  dynamic  study  is  also  carried  out,  three  additional 
scans  would  be  done: 

4.  Pre-contrast  3D  SPGR  scan:  the  proton  density-weighted  scan  is  done 


before  the  injection  of  the  contrast  material.  The  resulting  images  would  be 
used  to  compare  with  resulting  images  after  contrast  material  has  reached  a 
steady  state. 
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5.  2D  fast  SPGR  multi-slice  dynamic  contrast  study:  fourteen  identical  25- 
second  scans  are  performed  continuously  with  an  inter-scan  time  of  5  sec. 
At  each  scan,  about  1 1  slice  locations  are  imaged.  Three  scans  are 
performed  before  the  injection  of  the  contrast  material.  The  contrast- 
material  injection  and  the  fourth  scan  start  at  the  same  time.  The  total  scan 
time  for  the  whole  dynamic  contrast  study  is  about  7  minutes. 

6.  Post-contrast  3D  SPGR  scan:  the  same  scan  as  Scan  4  is  done  after  the 
injection  of  the  contrast  material.  The  effect  of  the  contrast  material  has 
reached  a  steady  state.  The  resulting  images  would  be  used  to  compare 
with  those  obtained  before  the  contrast  material  was  administrated. 

Examples  of  the  scan  procedure  and  the  details  of  the  scan  protocols  used  for  the 
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breast  studies  are  shown  in  Appendix  C.  In  the  following  section,  three  case  studies 
with  normal  breasts,  breasts  with  benign  abnormal  masses,  and  breasts  with  malignant 
tumors  will  be  discussed. 

8.3.1  Normal  Breast  Study 

In  this  study,  no  first-pass  dynamic  contrast  enhanced  study  was  performed. 
Ten  slice  locations  were  studied.  Figure  8.6  shows  images  at  one  slice  location  at  the 
Selective,  Non-selective  and  Regular  conditions  in  the  arterial  spin  tagging  scan.  As 
shown  in  this  figure  as  well  as  in  clinical  images,  the  fat  and  non-fat  regions  are  well 
separated  geometrically.  The  non-fat  region  is  at  the  inner  part  of  the  breasts  while  the 
fat  region  is  at  the  outer  part.  Based  on  the  nice  separation,  regions  of  interest  (ROIs) 
can  be  easily  isolated  for  the  characteristic  comparison  between  fat  and  non-fat  areas 
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:um  at  151. 


of  the  breast.  The  ROIs  with  the  size  of  35  mm  have  been  identified  from  all  10  slices 
acquired  from  the  arterial  spin  tagging  sequence.  Figure  8.7  shows  the / IX  and  T\n 
distribution  for  this  study,  and  can  be  summarized  in  Table  8.2  below: 

Table  8.2.  The/ IX  and  T\n  Summary  for  a  Study  with  Normal  Breasts 


Mean  of  Tin 
(sec) 

STD  of  Ti„ 
(sec) 

Mean  of 
f/X  (sec'1) 

STD  of 
/  IX  (sec'1) 

Fat 

0.435 

0.207 

-  0.005 

0.135 

Non-fat 

1.446 

0.403 

0.033 

0.044 

8.3.2  Breast  Study  with  a  Benign  Abnormal  Mass 

No  first  pass  contrast  enhanced  study  was  performed  in  this  study  either.  Ten 
slice  locations  were  studied.  Figure  8.8  shows  images  at  one  slice  location  from 
clinical  scans  as  well  as  from  the  arterial  spin  tagging  scan.  A  rod-shape  abnormal 
mass  was  easily  identified  in  the  clinical  scans  (Figure  8.8  (d)  and  (e)).  This  mass  was 
considered  benign  because  the  volunteer  aclcnowledged  that  it  had  been  there  and  had 
remained  the  same  size  for  more  than  10  years.  ROIs  with  the  size  of  35  mm2  have 
been  isolated  from  the  edge  and  the  center  of  the  abnormal  mass  as  well  as  non-fat 
area  outside  of  the  abnormal  mass  for  the  MR  tissue  characteristic  comparison.  Figure 
8.9  shows  the  f/X  and  T\„  distribution  for  these  different  tissues,  and  can  be 
summarized  in  the  table  following  (Table  8.3): 
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fIX  and  T’ln  Distribution  for  Normal  Breast  Tissue 


Figure  8.7.  The / IX  and  Tln  distribution  for  one  of  the  studies  with 
normal  breasts. 
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displayed  at  the  same  widow  leve. 
with  the  maximum  at  390  and  the 
minimum  at  195. 
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/ IX  and  rin  Distribution  for  Benign  Breast  Lesion 
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Figure  8.9.  The/M  and  r,n  distribution  for  one  breast  study 
with  a  benign  lesion. 
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Table  8.3.  The / IX  and  Tin  Summary  for  a  Breast  Study 
with  a  Benign  Abnormal  Mass 


Mean  of 

Pm  (sec) 

STD  of  Tin 
(sec) 

Mean  of 
fIX  ( sec'1) 

STD  of 
flX{  sec'1) 

Lesion  Center 

0.813 

0.053 

0.057 

0.018 

Lesion  Edge 

0.615 

0.140 

0.034 

0.018 

Non-fat  Tissue 

1.142 

0.487 

0.013 

0.025 

8.3.3  Breast  Study  with  a  Malignant  Tumor 

No  first  pass  contrast  enhanced  study  was  performed  in  this  study  either.  Three 
slice  locations  centering  around  the  abnormal  mass  were  studied.  The  malignant 
regions  were  identified  based  on  the  image  processing  and  visualization  techniques 
discussed  in  Chapter  7  and  were  labeled  by  color  mapping  onto  the  clinical  images 
(Figure  8.10).  These  regions  were  close  to  the  abnormal  mass  that  could  be  palpated, 

i 

and  its  malignancy  was  later  confirmed  by  needle  biopsy.  ROIs  with  the  size  of  35 
mm  has  been  isolated  at  the  tumor,  non-fat  regions  outside  the  tumor  as  well  as  fat 

tissue  for  MR  tissue  characteristic  comparison.  Figure  8.1 1  shows  the / IX  and  T\n 

\ 

distribution  for  these  different  areas,  and  can  be  summarized  in  the  table  below  (Table 
8.4): 


Table  8.4.  The  fIX  and  T\n  Summary  for  a  Breast  Study 
with  a  Malignant  Tumor 


Mean  of  Tin 
(sec) 

STD  of  Tm 
(sec) 

Mean  of 
fIX  (sec'1) 

STD  of 
/// l  ( sec'1) 

Fat 

0.385 

0.085 

-0.177 

0.592 

Non-fat 

0.842 

0.392 

0.035 

0.266 

Tumor 

1.145 

0.372  . 

0.111 

0.142 

As  shown  in  Figure  8.1 1,  because  fat  has  Tin's  much  lower  than  those  of  the  other  two 
groups  of  tissue,  it  can  be  separated  out  easily  based  on  Tin  value  alone.  The  mean 
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Proton  Density  Weighted,  Slice  6  T2  Weighted,  Slice  6 


Proton  Density  Weighted,  Slice  7  T2  Weighted,  Slice  7  , 


Proton  Density  Weighted,  Slice  8 


T2  Weighted,  Slice  8 


Figure  8.10.  Suspicious  cancer  region  color  mapping.  The  identified 
suspicious  cancer  pixels  are  mapped  onto  the  high-resolution  traditional 
clinical  MR  images  according  to  their  corresponding  slice  locations.  These 
images  are  from  three  consecutive  3mm  slices  with  inter-slice  distance  of 
10  mm.  The  white  mapping  pixels  are  actually  red  in  the  visualization 
program  BreastView. 
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/ IX  and  Tln  Distribution  for  Breast  Tissue  with  a  Malignant  Tumor 


Tin  (sec) 


A:  Fat  '  +:  Non-fat  o:  Tumor 


(a) 

Figure  8.1 1.  The  f  IX  and  Tln  distribution  for  breast  tissue  with  a  malignant 
tumor.  Figure  (b)  is  a  magnification  of  Figure  (a). 
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f IX  and  T[n  Distribution  for  Breast  Tissue  with  a  Malignant  Tumor 
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Figure  8.11.  (continued) 
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and  standard  deviation  of  the / IX  for  fat  have  a  wide  range,  and  thus / IX  does  not 
provide  any  value  for  fat  tissue  recognition.  Since  tumor  tissue  is  non-fat  tissue  and 
should  have  7jn  value  at  the  range  of  healthy  non-fat  tissue,  the  overlap  between  fat 
and  tumor  is  minimal.  However,  to  separate  tumor  tissue  from  non-fat  tissue,  both  the 
T\n  and / IX  values  should  be  used.  For  example,  as  shown  in  Table  8.4  and  Figure 
8.11,  both  the  7jn  and / IX  values  at  tumor  tissue  are  relatively  higher  than  those  at 
non-fat  tissue,  but  with  some  overlap. 

8.3.4.  Patient  Study  Summary 

Eighteen  subjects,  including  1 1  patients  referred  by  their  physicians  and  seven 
volunteers  referred  by  themselves,  participated  in  the  breast  studies.  Three  patients  and 
two  volunteers  were  eliminated  from  the  analysis  because  either  the  studies  were  not 
done  according  to  protocol,  there  was  excessive  motion  during  scanning,  the  abnormal 

V 

mass  was  too  close  to  the  chest  wall  to  be  measured  reliably,  or  no  patient  biopsy  was 

v. 

performed  to  confirm  the  results  of  the  MR  study.  There  were  a  total  of  13  useful 

cases  left  for  analysis.  For  all  these  cases,  the  suspicious  pixels  were  identified  based 

on  the  following  criteria  as  was  suggested  in  Chapter  7: 

7jn:  >  0.5  sec 
f!X\  >0.1  sec'1 
STD  of f/X\  <  0.1  sec'1 

and  the  suspicion  level  threshold  was  set  at  20.2%.  Analysis  has  been  done  using  the 
program  BreastView.  Appendix  D  shows  three  case  studies  using  this  program.  For 
information  on  how  to  use  this  program,  refer  to  Appendixes  B  and  G.  Table  8.5 
shows  the  result  of  these  13  cases. 
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Table  8.5.  Summary  of  Breast  Case  Studies 


Date  of 
Study 

Subject 

Referred 

by 

Palpable 
(Yes  or 
No)? 

Seen  on 
Clinical 
MRI  ? 

Level  of 
Suspicion 
Based  on 

T\n  and  f/X 

Assessment 
Based  on 
Suspicion 
Level 

Abnormal  Mass 
Type 

05/03/97 

Physician 

No 

Yes 

Low 

Negative 

Benign 

(Fibroadenoma) 

05/23/97 

Physician 

Yes 

Indeter¬ 

minable 

Moderate 

High 

Positive 

Malignant 

(Ductal 

Carcinoma) 

08/01/97 

Physician 

Yes 

No 

High 

Positive 

Malignant 

(Ductal 

Carcinoma) 

11/15/97 

Physician 

Yes 

No 

Low 

Negative 

Benign 

12/10/97 

Physician 

Yes 

No 

High 

Positive 

Benign 

04/13/98 

Physician 

Yes 

Indeter¬ 

minable 

Moderate 

High 

Positive 

Benign 

07/27/98 

Physician 

Yes 

No 

Low 

Negative 

Benign 

08/07/98 

Physician 

Yes 

No 

High 

Positive 

Benign 

02/07/97 

Herself 

No 

No 

Low 

Negative 

Normal 

03/05/97 

Herself 

Yes 

Yes 

Low 

Negative 

Benign 

(Fibroadenoma) 

07/26/97 

Herself 

No 

No 

Low 

Negative 

Normal 

11/16/97 

Herself 

Yes 

Yes 

Low 

Negative 

Benign 

(Fibroadenoma) 

12/14/97 

Herself 

No 

No 

Low 

Negative 

Normal 

Normal  =  normal  subject  who  does  not  have  any  tumor  based  on  reliable  health 
history, 

Abnormal  Mass  Type  =  determined  based  on  biopsy  or  reliable  health  history. 

Based  on  the  summary  shown  in  Table  8.5,  the  number  of  true  positive  (TP)  is 
2;  the  number  of  false  negative  (FN)  is  0;  the  number  of  false  positive  (FP)  is  3  and 
the  number  of  true  negative  (TN)  is  8.  Therefore,  the  true-positive  fraction  or 
sensitivity  is  100%,  calculated  by  TP/(TP+FN).  The  false-positive  fraction  is  27.3%, 
calculated  as  FP/(FP+TN).  The  specificity  is  72.7%,  which  is  equal  to  one  minus  the 
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false-positive  fraction.  Because  the  number  of  cases  investigated  was  limited,  these 
statistical  figures  cannot  provide  enough  confidence  in  evaluating  the  technique.  If 
more  breast  cases  have  been  studied,  a  ROC  curve,  commonly  used  to  evaluate  the 
performance  of  a  technique  in  radiology,  can  be  built  by  using  different  suspicion 
level  thresholds  or  different  thresholds  of  T\n,  f IX  and  STD  of  f  IX.  If  a  large  pool  of 
case  studies  is  available,  a  more  sophisticated  Bayesian  statistical  model  discussed  in 
Chapter  7  can  also  be  built.  This  model  can  be  used  to  evaluate  the  imaging  technique, 
and  more  importantly,  it  can  be  used  to  evaluate  new  case  studies. 


Chapter  9 

Summary,  Conclusions  and  Recommendations 


The  arterial  spin  tagging  technique  based  on  2D  fast  SPGR  has  been  developed 
and  implemented  in  both  EPIC  5.4  for  the  Signa  Advantage  MR  system  and  EPIC  LX 
8.2.5  for  the  Signa  Horizon  LX  MR  system.  Its  accuracy  in  perfusion  measurement 
has  been  proven  through  phantom  experiments.  It  has  also  been  successfully  applied  to 
tissue  classification  and  tumor  identification.  Its  applicability  is  enhanced  through  the 
graphical-user-interface  analysis  software,  BreastView,  which  has  been  developed  by 
the  author  for  the  necessary  image  visualization,  threshold  analysis,  and  cross¬ 
registration  between  different  MR  breast  imaging  techniques.  Based  on  the  13  case 
studies  analyzed,  a  sensitivity  of  100%  and  a  specificity  of  72.7%  have  been  achieved 
using  the  fast  SPGR  based  arterial  spin  tagging  sequence  and  the  associated  analytical 
techniques.  s 

However,  the  fast  SPGR  based  spin  tagging  technique  as  well  as  the  data 
analysis  technique  used  has  limited  the  accuracy  of  the  T\  and / /A  estimation  due  to 
two  ideal  assumptions.  First,  it  has  been  assumed  that  data  acquisition  would  not 
disturb  the  recovery  curves  of  the  inverted  spins  or  the  equilibrium  states  of  the  non- 
inverted  spins.  This  assumption  indeed  never  holds.  Simulation  showed  that  the  T\ 
calculated  based  on  the  techniques  presented  was  slightly  lower  than  the  true  T\.  The 
better  approach  is  to  replace  the  2D  fast  SPGR  data  acquisition  scheme  with  an  echo 
planar  data  acquisition  scheme.  In  echo  planar  imaging  (EPI),  only  one  read-out  RF 
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pulse  is  used  for  the  complete  data  acquisition  of  a  whole  image.  The  disturbance  due 
to  the  read-out  RF  pulses  to  the  recovery  curves  of  the  inverted  spins  or  the 
equilibrium  states  of  the  non-inverted  spins  would  be  small.  Since  EPI  is  faster  than 
2D  fast  SPGR,  the  total  scan  time  would  be  shortened  also.  This  would  enhance  its 
value  as  a  breast-screening  tool.  However,  EPI  places  high  demands  on  gradient 
amplitudes  and  slew  rates,  and  introduces  hard-to-remove  N/2  ghost  artifacts.  Thus  the 
2D  fast  SPGR  based  spin  tagging  sequence  can  still  play  a  useful  role  in  research  sites 
without  sufficient  hardware  capability.  This  EPI  based  spin  tagging  sequence  has  been 
implemented  successfully  in  the  simulation  level  in  EPIC  LX  8.2.5  for  the  Signal 
Horizon  LX  MR  systems.  It  has  not  been  able  to  run  successfully  on  the  LX  system 
yet  due  to  hardware  problems  that  GE  Medical  System  Inc.  is  still  attempting  to 
resolve. 

V 

The  other  assumption  for  the  2D  fast  SPGR  based  spin  tagging  sequence  is  that 
the  contribution  of  magnetization  transfer  in  Eq.  (4.1)  is  negligible.  RF  and  gradient 
pulses  dedicated  to  magnetization  transfer  have  been  implemented  in  the  EPI  based 
arterial  spin  tagging  pulse  sequence.  This  would  potentially  improve  the  calculation  of 
the  T\  values.  Relatively  higher  magnetization  transfers  have  also  been  detected  in 
malignant  tumor  [36].  Thus,  the  magnetization  transfer  parameters  detected  might  add 
another  feature  in  determining  lesion  malignancy.  Analysis  of  the  magnetization 
parameters  has  not  been  done  in  this  dissertation,  but  remains  as  a  recommendation  for 
the  further  development  of  this  project. 
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The  pulse  sequence  and  post-processing  techniques  presented  may  be  suitable 
for  all  specific  clinical  indications  that  have  been  reported  for  first-pass  contrast 
enhanced  imaging.  They  may  in  fact  be  more  convenient  and  equally  reliable  for 
identifying  the  regions  of  tissue  changes  that  give  rise  to  contrast  enhancement. 
However,  contrast  enhancement  is  regarded  as  essential  in  MRI  breast  evaluation,  and 
is  unlikely  to  be  supplanted.  The  recommendation  is  to  perform  the  first-pass  contrast 
enhancement  comparison  study  along  with  all  the  arterial  spin  tagging  studies.  This 
would  allow  the  direct  comparison  in  accuracy,  precision  as  well  as  ease  of 
application. 

For  the  further  development  of  this  project,  more  case  studies  using  all  the 
sequences  developed  should  be  conducted  to  develop  a  meaningful  statistical  model 
appropriate  for  Bayesian  analysis  as  discussed  in  Chapter  7.  To  make  it  successful, 
strong  collaboration  with  referring  physicians  and  surgeons  is  recommended.  This 
would  allow  sufficient  numbers  of  patient  volunteers.  Needle  biopsies  should  be 
performed  on  all  volunteers  with  palpable  masses  or  suspicious  masses  identified  in 
any  imaging  modality,  so  that  results  from  arterial  spin  tagging  studies  as  well  as  first- 
pass  contrast  enhanced  studies  can  be  confirmed. 
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Appendix  A 

Computer  Simulation  and  Scanner  Operation  for  the  Fast  SPGR  Based 
Arterial  Spin  Tagging  Pulse  Sequence 

The  fast  SPGR  based  arterial  spin  tagging  pulse  sequence  (BreFgre)  has  been 
developed  in  the  EPIC  LX  8.2.5  environment  on  an  SGI  workstation  with  an  operation 
system  of  IRIX  6.3  or  higher.  It  can  be  simulated  on  an  SGI  workstation  or  run  on  a 
Signa  Horizon  LX  system. 

A.l.  Simulation 

At  the  directory  containing  the  BreFgre  pulse  sequence  source  code,  run  the 
EPIC  tool  by  typing: 

t 

WTools.x  & 

Build  the  PSD  first  if  it  has  not  been  done  yet.  To  see  the  simulation,  click  the  “Load 

cv”  button  at  the  Evaltool  page,  and -then  load  the  file,  “GoodSimCVsl.txt”  or 

* 

“GoodSimCVs2.txt”.  The  file  “GoodSimCVsl.txt”  is  for  a  full  simulation.  The  file 

“GoodSimCVs2.txt”  is  for  a  fast  simulation  without  going  through  all  the  repetition 

steps.  Be  sure  to  press  the  return  key  after  each  entry  to  ensure  that  the  value  is 

accepted.  Follow  the  steps  for  simulation  at  the  IPGSim  page  to  view  the  waveform. 

The  CVs  from  these  two  files  only  serve  as  example  CVs;  they  can  be  changed  as 

desired.  The  CVs  used  in  the  two  simulation  files  are  defined  as  following: 

opfast  =  the  fast  scan  flag  with  1  =  on, 

BreScan  =  the  breast  scan  flag  with  1  =  on, 

opacqo  =  the  acquisition  mode  with  0  =  interleaved  and  1  =  sequential, 
opmph  =  the  flag  for  multi-phase  data  acquisition  with  1  =  on, 
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opfphases  =  the  total  images  acquired  at  each  slice  location  in  the  multi-phase 
case, 

opdeprep  =  the  flag  for  DE  prep  pulse  with  0  =  off  and  1  =  on, 
opirprep  =  the  flag  for  IR  prep  pulse  with  0  =  off  and  1  =  on, 
opirmode  =  the  type  of  IR  or  GRASS/MPGR  used, 
slquantl  =  the  total  number  of  images  acquired, 
phorder  =  the  phase/view  ordering  with  0  =  normal,  1  =  centric  and  2  = 
interleaved, 

br  norm  =  the  number  of  images  acquired  at  the  regular  condition, 
br  tino  =  the  number  of  TI  choices  used, 
ti_inc  =  the  step  size  of  TI  increment  (ms), 

br  acqs  =  the  number  of  images  repeatedly  acquired  at  each  TI  value  at  the 
spin  tagging  condition. 

A.2.  Scanner  Operation 

The  breast  sequence  scans  only  at  one  slice  location  at  a  time.  After  the  scan  is 
finished  at  each  location,  the  series  should  be  copied,  and  the  slice  location  should  be 
updated.  After  this,  click  through  the  sequence  of  “Save  Series”,  “Prepare  to  Scan” 

t 

and  "Scan".  It  is  not  necessary  to  do  the  full  prescan  at  every  slice  location;  every 
other  slice  location  is  sufficient.  This  reduces  the  overall  scan  time.  To  skip  the 
prescan,  run  "manual  prescan",  but  hbort  it  prematurely,  so  that  the  old  prescan  values 
will  be  kept. 

The  following  is  a  complete  example  prescription  for  a  breast  study  using  the 

BreFgre  sequence: 

Patient  ID:  0000027 
Patient  Name:  amy  doe 
Landmark:  Xyphoid 
Patient  position:  Prone 
Patient  Entry:  Head  First 
Coil:  BREASTPA 

Series  Description:  david's  breast  seq 
Plan:  coronal 
Mode:  2D 
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Pulse  Seq:  SPGR 
Imaging  Options:  Seq,  Fast,  Mph 
Psd  Name:  /psd/dzhu/BreFgre 
TE:  Minimum 
Flip  Angle:  10° 

Bandwidth:  31  kHz 

FOV:  34  cm 

Slice  Thickness:  3.0  mm 

Spacing:  1.0  mm 

Freq:  256 

Phase:  256 

NEX:  1 

Phase  FOV:  0.5 
Freq  Dir:  R/L 
Auto  Center  Freq:  Water 
Start:  A44 
End:  A44 
#  Slices:  1 

At  Multi  Phase  Screen: 

Phase  Per  Location:  61 
Phase  Acq.  Order:  Sequential 
Delay  After  Acq.:  Minimum  1 

After  all  the  necessary  CVs  have  been  entered,  click  through  the  sequence  of 
“Save  Series”,  “Prepare  to  Scan”  and  "Scan".  After  the  scan  is  finished,  the  images 
would  be  reconstructed  and  displayed  on  the  console. 


Appendix  B 

Instructions  for  the  BreastView  Program 


The  BreastView  program  has  been  designed  to  be  user-friendly.  The  graphical- 
user-interface  guides  the  user  through  the  entire  visualization  and  analysis  process  for 
a  breast  study.  The  assumption  is  that  the  necessaiy  analysis  has  been  done  correctly 
beforehand.  After  executing  the  BreastView  program  on  the  command  line  by  typing 
"BreastView",  there  is  a  "Help"  button  on  the  first  page  of  the  interface.  Clicking  the 
"Help"  button  allows  the  user  to  view  the  full  html  file  with  the  following  instructions: 

Installation 

At  the  program  directory  of  BreastView,  create  the  following  symbolic  linkage: 

In  -s  /CDROM  DatalnCD 

In  -s  /data2/zhudata/BreastData  DatalnOther 

In  -s  /data2/zhudata/BreastAnaly  AnalyDir 

To  run  BreastView,  type 
BreastView  & 

Data  Loading 

Click  the  "Load  Data  Directory"  button  to  open  the  directory  where  data  is  stored.  The 
data  is  in 

DatalnCD 

or 

DatalnOther 
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After  the  data  directory  is  loaded,  the  "Show  Study  Info"  and  "Load  Analysis 
Directory"  buttons  will  come  up.  Clicking  the  "Show  Study  Info"  button  will  show  the 
information  of  the  breast  study  that  has  been  done. 

Click  the  "Load  Analysis  Directory"  button  to  open  the  directory  where  analytical 
results  are  stored.  The  analysis  directory  should  be  in 

AnalyDir 

If  this  directory  has  been  loaded  correctly,  all  the  action  buttons  will  come  up. 
Otherwise  the  program  will  ask  you  to  reload  the  analysis  or  data  directory. 

Data  Preview 

Clicking  the  "Display  Histogram  for  All"  button  displays  the  histogram  of  T\  and / //L 
for  all  the  slices  investigated. 

Create  Suspicious  Cancer  Regions 

\ 

The  detection  images  have  normally  been  created  according  to  the  default  threshold 
settings.  They  can  be  recreated  based  on  the  researcher's  choice  of  thresholds.  After 
the  detection  thresholds  have  been  specified,  simply  click  the  "Create  Detection 
Image"  button  to  create  a  new  set  of  detection  images. 

View  All  Images 

Clicking  the  "Display  All  Images"  button  should  lead  to  the  display  of  the  eight  types 
of  images  that  are  most  useful  for  making  decisions. 

This  button  offers  two  choices: 
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1)  Contrast  enhancement  dynamic  study  images  will  be  shown. 

2)  No  contrast  enhancement  dynamic  study  images  will  be  shown.  All  feature 

images  will  be  displayed  instead. 

Tips  in  Viewing  Images 

After  all  eight  types  of  images  are  displayed,  the  user  can  do  the  following: 

1)  The  Maximum  and  Minimum  values  of  each  image  are  displayed  automatically. 

The  image  contrast  can  be  changed  by  moving  the  "Max"  and  "Min"  bars. 

2)  Slices  of  investigation  can  be  displayed  by  moving  the  sliding  bars.  By  moving 
the  sliding  bars,  the  slice  locations  of  different  types  of  study  can  be  approximately 
matched. 

3Dslice#:  Select  the  slice  of  the  3D  images. 

DySlice#:  Select  the  slice  for  the  dynamic  study,  and  the  corresponding 
slices  of  the  3D  studies. 

BrSlice#:  Select  the  slice  for  all  the  feature  images  resulted  from  the  spin 
tagging  sequence  and  the  slices  of  the  corresponding  dynamic  studies. 

HiSlice#:  Select  the  slice  of  the  high-resolution  clinical  images  (Proton- 
weighted  and  72-weighted). 

AxSlice#:  Select  the  slice  of  the  axial  view  images. 

3)  To  determine  an  individual  pixel  value,  simply  click  the  left  mouse  button  on 
the  pixel  of  interest  on  any  image  except  the  axial  image.  Other  images  should 
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display  the  pixel  values  at  the  corresponding  pixel  location  also.  The  units  of  the 
pixel  values  are  indicated  if  known. 

4)  Click  the  "Suspicion  mapping"  button  to  show  the  color  mapping  in  clinical  and 
dynamic  rise  images.  The  color  mapping  corresponds  with  the  range  setting  for  the 
Suspicion  map.  Changing  the  range  setting  of  the  Suspicion  map  through  the 
"Max"  and  "Min"  sliding  bars  should  automatically  display  the  new  color 
mappings. 

5)  Click  the  "Plot"  button  to  open  a  window  for  plotting  the  dynamic  study. 

Select  "Point"  for  pixel  by  pixel  display.  Clicking  the  left  mouse  button  on  a  pixel 

in  any  image  except  the  axial  image  will  show  the  dynamic  plot. 

Select  "Roi"  to  show  the  dynamic  plot  of  the  ROI  (region  of  interest)  specified 
earlier. 

6)  After  activating  "Draw  ROI ...",  the  ROI  can  be  specified  by  the  user:  left 
mouse  button  to  specify  the  starting  point  of  the  ROI,  right  mouse  button  to  create 
the  ROI,  and  middle  button  to  finish  the  task. 

7)  If  the  dynamic  study  images  are  displayed,  you  have  a  choice  of  choosing  the 
"Axial  Cut"  or  the  "Dynamic  Rise"  images. 

8)  Click  the  "Histogram"  button  to  display  the  total  histogram  including  all  the 
slices  investigated. 


Pg .138 


Click  the  "Pixel"  button:  the  T\  and / IX  values  will  be  mapped  onto  the  total 
histogram  while  clicking  the  left  mouse  button  on  the  images  displayed. 

Click  the  "Roi"  button  to  map  the  distribution  of  7j  and / IX  of  the  ROI  specified 
on  the  total  histogram. 

Click  the  "No  Mapping"  button  to  stop  the  above  mapping. 

9)  The  slice  selection  bars  have  been  designed  to  view  the  appropriate  images.  By 
dragging  them,  the  user  can  see  how  they  work. 

10)  Click  the  "Redraw  All  Images"  button  to  redraw  all  the  images. 

1 1)  An  active  window  can  be  closed  by  selecting  "Close"  from  the  pull-down 

i 

manual  at  the  upper-left  comer. 

Close  Current  Study  and  Open  Another  One 

To  close  the  current  study,  just  click  the  "Close  Study"  button,  then  you  are  ready 
to  analyze  another  study. 
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Appendix  C 

Breast  Study  Procedures  and  Protocols 


Two  example  procedures  with  and  without  a  first-pass  contrast  enhanced 
dynamic  study  are  shown  in  Sections  C.l  and  C.2.  The  investigator  must  write  down 
the  relevant  information  during  the  study.  If  there  is  an  abnormal  mass  for  the  subject, 
the  subject  or  her  physician  should  be  asked  to  localize  the  masses  on  the  figure 
provided  (as  shown  in  Section  C.3).  The  other  protocols  are  shown  in  Sections  C.4  to 
C.8. 

C.l.  Example  Procedure  without  a  First-Pass  Contrast  Study 


Date  of  study:  1 1/16/97 

Breast  Sequence  Study 

by  Michael  Buonocore  and  David  Zfiu 

Type  of  study:  phantom( ),  subject  ( x),  Patient  ( ) 

Name:  Age:  38 

Preparation: 

Ask  subject  to  have  a  good  fit  of  the  foam  pads  and  a  good  strap,  not  to  move,  and 
breath  smoothly  during  scanning. 

Procedures: 

1 .  2D  Spin  Echo  axial  cut  (Series  1)  FOV  =  40  cm 
I  84  to  s  84  ,  loc:  29  .  Scan  thick  =  5  mm,  Inter  =  1  mm. 
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2.  2D  Fast  Spin  Echo  (Series  2)  FOV  =  34  cm 

P  7  to  A  55  ,  loc:  12  ,  scan  thick  =  5  mm,  inter  =  1  mm. 

3.  David's  breast  sequence  coronal  cut:  fgtl_new3d,  FOV  =  34  cm 
Seven  TI  cycles:  br_TR  =  2.7  s,  br_TI  =15  ms,  ti_inc  =100  ms. 

TI  decreases  from  615  ms  to  15  ms. 

Acquisition  slice  thickness:  3  mm,  Selective  slice  thickness:  3mm 
Frequency  direction:  R/L  for  all  locations 
series  3:  location:P7  R1  =  6  ,  R2  =  15  ,  TG  =153 

series  4:  location:Pl  R1  =  6  ,  R2  =  15  ,  TG  =154 

series  5:  location: A5  R1  =  6  ,  R2  =  15  ,  TG  =  153 

series  6:  location: A1 1  R1  =  6  ,  R2  =  1 5  ,  TG  =  1 54 

series  7:  location: A1 7  R1  =  6  ,  R2  =  15  ,  TG  =  160 

series  8:  location: A23  R1  =  6  ,  R2  =  15  ,  TG  =  160 

series  9:  location: A29 v  R1  =  6  ,  R2  =  1 5  ,TG=160 

% 

series  10:  location:A35  R1  =  6  ,  R2  =  15  ,  TG  =  167 
Comments: 

Fibroadenoma  in  the  left  breast,  about  1/3  in  from  nipple,  and  about  3  o'clock 
position. 
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C.2.  Example  Procedure  with  a  First-Pass  Contrast  Study 


Date  of  study:  8/7/98 

Breast  Study  with  Contrast 

by  Michael  Buonocore  and  David  Zhu 

Type  of  study:  phantom( ),  normal  subject  ( ),  subject  w/  suspicious  lesion  (x) 

Name:  Age:  37  Weight:  135 

Preparation: 

Ask  subject  to  have  a  good  fit  of  the  foam  pads  and  a  good  strap,  not  to  move,  and 
breath  smoothly  during  scanning. 

Procedures:  (Protocol:  #  64)  (Study  #  04817) 

i 

1 .  2D  Spin  Echo  axial  cut  (Series  1  ->  axial)  FOV  =  40  cm 

I  84  to  s  84  ,  loc:  29  .  Scan  thick  =  5  mm,  Inter  =  1  mm. 

2.  2D  Fast  Spin  Echo  (Series  2  ->  highres)  FOV  =  34  cm 

P  10  to  A  50  ,  loc:l  1  ,  scan  thick  =  5  mm,  inter  =  1  mm. 

3.  David's  breast  sequence  coronal  cut:  fgtl_new3d,  FOV  =  34  cm 
Seven  TI  cycles:  br_TR  =  2.7  s,  br_TI  =  15  ms,  ti_inc  =100  ms. 

Tl  decreases  from  615  ms  to  15  ms. 

Acquisition  slice  thickness:  3  mm,  Selective  slice  thickness:  3  mm 
Frequency  direction:  R/L  for  all  locations 
series 3:  location:  P 10  R1=6,R2  =  15  ,TG=146 

series  4:  location:  P4  R1  =  ,  R2  =  ,  TG  = 
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series  5:  location:  A2  R1  =  6 ,  R2  =  15  ,  TG  =  150 

series  6:  location:  A8  R1  =  6  ,  R2  =  15  ,  TG  =  149 

series  7:  location:  A14  R1  =  ,  R2  =  ,  TG  = 

series  8:  location:  A20  R1  =  6  ,  R2  =  15  ,  TG  =  146 

series  9:  location:  A26  R1  =  ,  R2  =  ,  TG  = 

series  10:  location:  A32  R1  =  6  ,  R2  =  15  ,  TG  =  149 

series  11:  location:  A38  R1  =  ,  R2  =  ,  TG  = 

series  12:  location:  A44  R1  =  6  ,  R2  =  15  ,  TG  =  158 
series  13:  location:  A50  R1  =  6,R2  =  15  ,TG=151 

4.  Pre-contrast  3D  scan  (scan  #2  in  protocol,  2'37")  (flip  angle  =  65) 

(series  14  ->  pre_3d_a) 

(Center  at  A22.7  with  28  Iocs,  3  mm  slice  thickness) 

5.  More  David's  breast  sequence  coronal  cut:  fgtl_new3d,  FOV  =  34  cm 

series  15:  location:  A23  R1  =  6  ,  R2  =  15  ,  TG  =  145 

series  16:  location:  A17  R1  =  ,  R2  =  ,  TG  = 

series  17:  location:  A29  R1  =  6  ,  R2  =  15  ,  TG  =  148 

The  patient  is  pulled  out  to  set  up  contrast  injection  in  the  following  scans: 

6.  Pre-contrast  3D  scan  (scan  #2  in  protocol,  2'37")  (flip  angle  =  65) 

(series  18  ->  pre_3d) 

(Repeat  the  pre_contrast)  (Center  at  A22.7  with  28  Iocs,  3  mm  slice 

thickness) 

7.  Fast  SPGR  contrast  dynamic  study  (total  7  minutes)  (scan  #1  in  protocol): 
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(Series  19  ->  dynamic) 

Each  scan  has  a  scan  time  of  25  sec,  with  an  approximate  inter-scan  time  of  5 
sec,  for  a  total  of  30  sec. 

(1) 3  scans  are  done  before  injection. 

(2)  Begin  injection  and  press  SCAN  to  begin  the  4th  scan. 

(3)  Continue  on  the  same  fashion  until  the  14th  scan. 

8.  Post-contrast  3D  scan  (scan  #2  in  protocol,  2'37"),  same  as  #6  above. 

( Series  20  ->  post_3d) 

(Center  at  A22.7  with  28  Iocs,  3  mm  slice  thickness) 

Comments: 

Patient  claims  that  the  tumor  is  located  at  10  o'clock  in  the  right  breast. 

The  patient  has  small  breasts.  Foam  packs  have  been  adjusted  to  avoid  motion 
artifacts.  Might  need  to  consider  image  displacement  correction. 

Based  on  the  analysis  of  the  perfusion  study,  this  patient  appears  to  have  a 
malignant  tumor  in  the  right  breast. 

Biopsy  later  indicated  that  the  lesion(s)  was  benign. 
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C.3.  A  Simple  Questionnaire  for  tbe  Localization  of  the  Abnormal  Masses 


Volunteer  name: 


Date  of  study: 


View  from  side 


C.4. 2D  Spin  Echo  Scan  in  the  Axial  Plane 

Patient  Entry:  head  first 
Patient  Position:  prone 
Landmark:  sternal  notch 
Coil:  other:  LS45 
Plane:  axial 

Mode:  2D  Monitor  SAR 
Pulse  Sequence:  spin  echo 
Resp  Flow 
comp  comp 
Number  of  Echo  =  1 
TE  =  Minimum  Full  (20  ms) 

TR  =  600  msec 
Autoshim 

Auto  Center  Freq:  Water 
FOV  =  40  cm 
Slice  Thickness  =  5.0  mm 
Inter-slice  Gap  =  1  mm 

Start:  184,  End:  S84  =>  Number  of  Locations:  29 
Acq  Time: 

Freq:  256 
Phase:  128 
Freq  Dir:  R/L 
Phase  FOV:  40  cm 
1  NEX:  3  min  38  sec 

C.5. 2D  Fast  Spin  Echo  in  the  Coronal  Plane 

Patient  Entry:  head  first 
Patient  Position:  prone 
Landmark:  sternal  notch 
Coil:  other:  LS45 
Plane:  Coronal 
Mode:  2D  Monitor  SAR 
Pulse  Sequence:  Spin  Echo 
Fast 

FSE  Optimization:  1 
Echo  Train:  8 
Number  of  Echo:  2 
EffTE:  17  ms 
Eff  TE2:  128  ms 
Rep  Time:  3600  ms 


Autoshim  Off  (only  need  to  do  once  for  the  whole  study) 

Auto  Center  Freq:  Water 
FAT  SAT 
FOV  =  34  cm 
Slice  Thickness  =  6  mm 
Inter-slice  Gap  =  1  mm 
Suggest  Range: 

Start:  P10,  End:  A62  =>  Number  of  Locations:  13 
Acq  Time: 

Freq:  256 
Phase:  192 
Freq  Dir:  R/L 
1  NEX:  3  min  12  sec 
Reps  before  pause:  None 

C.6.  Fast  SPGR  Based  Arterial  Spin  Tagging  Pulse  Sequence  in  the  Coronal 

Plane  in  Signa  Advantage  System 

Patient  Entry:  head  first 
Patient  Position:  prone 
Landmark:  sternal  notch! 

Coil:  other:  LS45 
Plane:  Coronal 

Mode:  2D  Turn  off  monitor  SAR 
Pulse  Sequence:  SPGR 
sequential 
fast  1 

multiphase 

Psd  Name:  zhud/fgtl_new3d 
Phase  Per  Location:  61 
Phase  Acq.  Order:  sequential 
Delay  After  Acq.:  minimum 
Flip  Angle:  10  degree 
Minimum  TE 
Autoshim  Off 
Auto  Center  Freq:  Water 
Bandwidth:  32  kHZ 
FOV  =  34  cm 
Slice  Thickness:  3.0  mm 
Inter-slice  Gap:  2.5  mm 

Start:  A30,  End:  A30  (or  other  locations)  =>  Number  of  Locations:  1 
Acq  Time: 
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Freq:  256 
Phase:  256 
Freq  Dir:  R/L 
Phase  FOV:  34  cm 
1  NEX 

Use  the  following  default  parameters  : 

br_tino  =  7,  br_TI  =  15ms,  br_TR  =  2.7  s,  br_opyres  =  120,  br_opyresh  =  120, 
yhybrid  =  1,  opfphases  =  61,  br_tino  =  7,  br_reps  =  1,  br_norm  =  5,  br_acqs  =  5, 
br_TI  =15  ms,  ti_inc  =100  ms. 

The  protocol  for  the  fast  SPGR  based  arterial  spin  tagging  sequence  in  Signa 

Horizon  LX  can  be  referred  to  Appendix  A. 

C.7. 3D  SPGR  Scan 

Patient  Entry:  head  first 
Patient  Position:  prone 
Landmark:  sternal  notch 
Coil:  other:  LS45 
Plane:  Coronal 
Mode:  3D  Monitor  SAR 
Pulse  Seq:  SPGR 
Flow  Comp 
Graphic  RX 
Flip  Angle:  io  degree 
TE  =  Minimum 
TR  =  50  msec 
Autoshim  Off 
Auto  Center  Freq:  Water 
FAT  SAT 
FOV  =  34  cm 
Slice  Thickness:  3.0  mm 
Number  of  Locations:  28 
Acq  Time: 

Freq:  256 
Phase:  128 
Freq  Dir:  R/L 

Phase  FOV:  3/4  =>  25.5  cm 
1  NEX:  2  min  37  sec 
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C.8.  First-Pass  Contrast  Enhanced  Dynamic  Study  using  Multi-slice  Fast  SPGR 

Patient  Entry:  head  first 
Patient  Position:  prone 
Landmark:  sternal  notch 
Coil:  other:  LS45 
Plane:  Coronal 
Mode:  2D  Monitor  SAR 
Pulse  Seq:  SPGR 
Fast 

Flip  Angle:  30  degree 

TE=  4.2  msec 

TR=  11.7  msec 

Autoshim  Off 

Auto  Center  Freq:  Water 

FOV  =  34  cm 

Slice  Thickness:  5  mm 

Inter-slice  Gap:  2  mm 

Start:  ,  End:  =>  Number  of  Locations:  11 

Acq  Time: 

Freq:  256 
Phase:  128 
Freq  Dir:  R/L 

Phase  FOV:  3/4  =>  25.5  cm  ‘ 

2NEX:  25  sec 

\  w 
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Appendix  D 

Example  Breast  Case  Studies 


For  all  these  cases,  the  suspicious  pixels  were  identified  based  on  the  following 

criteria  as  was  suggested  in  Chapter  7: 

7in:  >  0.5  sec 
f!X\  >0.1  sec'1 
STD  of f  IX:  <  0.1  sec'1 

and  the  suspicion  level  threshold  was  set  at  20.2%.  Analysis  has  been  done  using  the 
program  BreastView. 

Case  Study  #1:  Palpable  mpss 

Date  of  study:  Nov  16, 1997 

Level  of  suspicion  based  on  7i„  and  perfusion:  Low 

Clinical  MR  images;  Positive  (lesion  visualized,  consistent  with  palpation) 

Assessment  based  on  spin  tagging:  Negative  for  malignancy 

Abnormal  mass  type:  Benign  based  on  reliable  health  history  (CONSISTENT  with 

spin  tagging  result) 

Discussion:  The  low  level  of  suspicion,  based  on  T\n  relaxation  time  and  perfusion,  at 
and  around  the  palpated  lesion  location  was  a  strong  indication  that  the  lesion  was 
benign.  Only  one  suspicious  pixel  was  found  at  that  location  at  the  20.2%  threshold.  It 
was  judged  that  the  breasts  did  not  contain  malignant  tissue.  This  case  is  shown  in 
Figures  A.l(a)-(d). 
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Case  Study  #2:  Palpable  mass 

Date  of  study:  Aug  6, 1997 

Level  of  suspicion  based  on  Ti„  and  perfusion:  High 

Clinical  MR  images:  Negative  (mass  not  clearly  visualized,  prior  cyst  visualized) 
Assessment  based  on  spin  tagging:  Positive  for  malignancy 
Biopsy  result:  Positive  (spin  tagging  CONSISTENT  with  biopsy) 

Discussion:  The  high  level  of  suspicion,  based  on  Tin  relaxation  time  and  perfusion,  at 
and  around  the  lesion  detected  by  palpation,  was  a  strong  indication  that  the  lesion  was 
malignant.  This  case  is  shown  in  Figure  A.2(a)-(d). 

Case  Study  #3:  Palpable  mass 
Date  of  study:  Aug7,  1998 

Level  of  suspicion  based  on  T\n  and  perfusion:  High 
Clinical  MR  images:  Positive  (lesion  visualized) 

Assessment  based  on  spin  tagging:  Positive  for  malignancy 

Assessment  based  on  first-pass  contrast  enhancement  study:  Positive  for  malignancy 

Biopsy  result:  Negative  (spin  tagging  NOT  CONSISTENT  with  biopsy) 

Discussion:  High  levels  of  suspicion,  based  on  T\n  relaxation  time  and  perfusion,  were 
found  in  several  breast  regions.  The  visible  lesion  location  on  clinical  MRI  was 
consistent  with  that  of  the  palpable  mass.  It  overlapped  but  was  considerably  smaller 
than  that  found  by  the  levels  of  suspicion.  First-pass  contrast  enhancement  study 
results  correlated  strongly  with  the  high  perfusion  regions  detected  by  spin  tagging. 
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Nevertheless,  biopsy  showed  the  mass  to  be  benign.  This  case  is  shown  in  Figure 
A.3(a)-(d). 
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Figure  A.  1(a).  Nov  16,  1997  breast  study.  Suspicion  level  map 
is  created  using  the  “Create  Detection  Image”  button  after 
setting  the  thresholds. 
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Figure  A.  1(d).  Nov  16,  1997  breast  study.  The  Tx  and///l  values  of  pixels 
representing  the  lesion  identified  by  proton  density  and  T2  images,  are 
superimposed  (by  pure  white  pixels)  on  the  histogram  of  the  Tx  and / IX 
values  from  the  entire  breast.  This  figure  shows  that  the  criteria  for  suspicious 
pixels  is  based  on  moderate  Tx  and  high f IX.  See  text  for  complete  criteria. 
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Figure  A.2(a).  Aug  6,  1997  breast  study.  Suspicion  level  map 
is  created  using  the  “Create  Detection  Image”  button  after 
setting  the  thresholds. 
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Figure  A.2(d).  Aug  6, 1997  breast  study.  The  Tl  and f IX  values  of  pixels 
representing  the  lesion  identified  by  proton  density  and  T2  images  in  Figure 
A.2(c),  are  superimposed  (by  pure  white  pixels)  on  the  histogram  of  the  Tx 
and / IX  values  from  the  entire  breast.  This  figure  shows  that  the  criteria  for 
suspicious  pixels  is  based  on  moderate  Ti  and  high  f  IX.  See  text  for  complete 
criteria. 
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Figure  A.3(a).  Aug  7,  1998  breast  study.  Suspicion  level  map 
is  created  using  the  “Create  Detection  Image”  button  after 
setting  the  thresholds. 
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Appendix  F 

Some  Simple  and  Useful  Display  and  Analysis  Programs 

This  appendix  introduces  some  simple  but  useful  generic  image  display  and 
analysis  programs  that  have  been  written  by  the  author  or  adopted  from  other  sources. 
They  should  be  beneficial  for  any  research  work  in  MRI.  These  executable  files  and 
source  codes  are  available  in  the  CDROM  (zhusoft)  if  more  details  are  needed.  All 
these  programs  should  run  successfully  on  an  SGI  workstation  with  an  operating 
system  of  IRIX  6.3  or  higher. 

F.l.  dispAlls 

This  program  has  been  written  by  the  author  for  the  observation  and  ROI 

i 

analysis  of  one  image  or  a  series  of  images  without  a  header  or  with  a  header  but  in 
GE's  I.###  format.  To  see  how  to  use  this  pfogram,  simply  type  it  at  the  command  line 

without  an  argument: 

* 

>  dispAlls 

Functions:  (1)  Display  one  or  a  sequence  of  images  with  or  without  a  header, 

(2)  ROI  analysis  and  (3)  series  plot. 

Usage:  dispAlls  xres  yres  display_dim  image_name  headerStatus  MaxImagNum 
image_name  =  full  name  for  only  1  image  or  partial  name  for  a  sequence 
headerStatus  =  1:  with  header,  or  0:  without  header 
MaxImagNum  default  to  1  when  omitted 
Example:  dispAlls  256  256  256  i  1  29 
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Example:  dispAlls  256  256  256  is  0  80 
Example:  dispAlls  256  256  256  i.003  1  1 
Example:  dispAlls  256  256  256  is.003  0  1 
Example:  dispAlls  256  256  256  i.003  1 

F.2.  mydisplay  and  dispGE: 

These  two  programs  have  been  modified  from  Gao's  display_new  program 
[63],  They  are  simple  display  programs  for  one  image  with  or  without  a  header  file.  To 
see  how  to  use  them,  type  the  command  without  an  argument  as  following: 

>  mydisplay 

Display  a  no-header  image 

Usage:  mydisplay  xres  yres  display_dim  imagefile 

>  mydispGE 

Display  a  GE I.###  image 

Usage:  mydispGE  xres  yres  display J_dim  imagefile 

F.3.  Raw  Data  Inspection 

In  pulse  sequence  development,  raw  data  is  often  obtained  from  MR  scanning. 
The  raw  data  might  need  to  be  inspected  for  the  header  size,  the  center  of  k  space,  and 
whether  data  is  collected  correctly  in  the  frequency  domain.  This  purpose  can  be 
achieved  through  two  approaches.  The  first  approach  is  to  use  the  "od"  command  in  a 
UNIX  workstation: 

To  view  a  P  file  on  an  SGI: 

od  -s  Pfile  +  size_in_bytes_from_beginning_of_file.  |more 
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For  example, 

od  -s  Pfile  +39936.  jmore 
To  view  a  P  file  on  a  SUN: 

od  -s  Pfile  +  size_in_bytes_fi-om_beginning_of_file.  |more 
For  example, 

od  -s  Pfile  +39936.  |more 
or 

od  -iv  size_in_bytes_from_beginning_of_file.  |more 
For  example, 

od  -iv  Pfile  +39936.  |more 

The  second  approach  is  to  use  some  simple  Matlab  routines.  This  approach 

f 

allows  the  graphical  view  of  the  raw  data  and  should  provide  an  easier  understanding 
of  the  raw  data.  A  routine  written  by  the  author  is  shown  below  to  serve  as  an 
example:  ■ 

%  This  routine  is  to  view  the  P  file  contents 
clear; 

pfile  =  input('The  P  file  to  check  is: ',  's'); 

rawSize  =  input('The  size  of  the  P  files  (in  Bytes)  is: '); 

fid  =  fopen(pfile,'r'); 

[B,  count]  =  fread(fid,(rawSize/2),  'short');  % 

fclose(fid); 

count 
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beginP=  input('Set  the  beginning  point  for  the  P  file  viewing  (in  short): '); 
view__size  =  input('The  view  size  of  data  points  at  a  time: '); 
for  k  =  1 : 1000  %  give  a  large  number  1000  to  loop  forever 

start  =  input('The  chunck  number  with  the  view  size  just  set  or  zero  to  stop: '); 
start_in_short_at  =  (start- l)*view_size*2  +  beginP 
if  start  =  0 

break; 

else 

for  m  =  1 :  viewsize 

realB(m)  =  B(start_in_short_at  +  2*m  -1); 
imagB(m)  =  B(start_in_short_at  +  2*m); 

i 

end 

subplot(2,l,l); 

plot(realB); 

* 

title('Real  Data  Points'); 

subplot(2,l,2); 

plot(imagB); 

title('Imaginary  Data  Points'); 
end 
end 
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Appendix  G 

The  Companion  CDROM  Zhusoft 


The  companion  CDROM  Zhusoft  or  its  contents  can  be  obtained  from  the 
author,  David  Zhu,  or  his  major  professor,  Dr.  Michael  Buonocore  in  the  Department 
of  Radiology  at  University  of  California  Davis,  through  necessary  research 
agreements.  This  CDROM  contains  the  source  code  developed  for  this  dissertation.  It 
also  contains  some  executable  programs  that  can  be  run  directly  from  the  CDROM  on 
an  SGI  workstation  running  IRIX  version  6.3  or  higher.  A  full  breast  study  and  images 
from  a  4D  human  cardiac  study  have  also  been  included  so  that  analysis  can  be 
demonstrated  for  the  BreastView  program  and  the  dispAlls  program  respectively.  The 

t 

source  code  contained  in  this  CDROM  can  serve  as  examples  and/or  building  blocks 
for  the  further  development  of  the  breast  and  4D  imaging  projects  as  well  as  the 

development  of  MR  pulse  sequences  and  graphical-user-interface  analysis  applications 

* 

in  general.  Figure  A.4  shows  the  structure  of  software  and  data  contained  in  the 
CDROM.  Each  directory  contains  a  'Teadmefirst.txt"  text  file  that  explains  the 
contents  of  the  directory  and  gives  necessary  instructions. 

Under  the  directory  "4DpcExample",  the  subdirectory  "jul29"  contains 
resulting  images  from  a  4D  study.  They  can  be  viewed  and  analyzed  using  the 
dispAlls  program. 

The  directory  "/Breast/breastview/"  contains  the  breast  image  analysis 
program,  BreastView,  and  its  source  code.  Study  aug7  serves  as  an  example  for  the 
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analysis.  The  directory  "/Breast/Im_process/"  contains  the  entire  image  processing 
source  code  for  the  arterial  spin  tagging  imaging. 


CDROM  root  directory 


jul29  breastview  Im_process  Breast  4Dpc  EPI  Kidney 

/  \ 

DataDir  AnalyDir 

I  I 

aug7  aug7 


Figure  A.4.  The  data  and  software  structure  of  the  companion 
CDROM  Zhusoft. 


The  directory  "PulseSeqsAndRecons"  contains  all  of  the  MR  pulse  sequences 
and  the  necessary  image  reconstruction  programs  developed  for  this  dissertation.  The 
subdirectory  "Breast"  contains  the  pulse  sequences  for  the  breast  project.  The 
subdirectory  "4Dpc"  contains  the  pulse  sequences  for  the  4D  project.  This 
subdirectory  also  contains  the  two  4D  image  reconstruction  programs,  4dReconFrLx 
and  4dReconFr2AcqLX,  and  their  C  code.  Although  EPI  itself  and  the  study  of  kidney 
are  not  parts  of  this  dissertation,  the  author  has  been  involved  in  research  in  these 
areas.  Pulse  sequences  and  the  necessary  image  reconstruction  programs  developed  by 
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the  author  for  the  odd-number  hybrid  EPI  technique  and  for  the  study  of  kidney  are 
also  included  in  the  CDROM  as  an  extra  resource  under  the  subdirectories  "EPI"  and 
"Kidney"  respectively. 

The  directory  "DispProgs"  contains  all  the  generic  display  programs  that  have 
been  discussed  in  Appendix  F,  including  dispAlls,  mydisplay,  mydispGE,  and  some 
Matlab  source  code  for  raw  data  inspection. 

The  directory  "Dissertation"  contains  the  entire  dissertation,  which  can  be 
viewed  with  Microsoft  Word  97  and  Microsoft  PowerPoint  97  or  higher  versions. 
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